(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 11.2' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    158216,       3288]
NotebookOptionsPosition[    148749,       3172]
NotebookOutlinePosition[    149113,       3188]
CellTagsIndexPosition[    149070,       3185]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Horan", " ", "and", " ", 
    RowBox[{"Shortle", ".", " ", "The"}], " ", "Song", " ", "Remains", " ", 
    "Not", " ", "the", " ", 
    RowBox[{"Same", ":", " ", 
     RowBox[{
     "Correlated", " ", "Intercept", " ", "and", " ", "Slope", " ", 
      "Uncertainties", " ", "Matter", " ", "to", " ", "Prices", " ", "vs", 
      " ", 
      RowBox[{"Quantities", ".", "  ", "JAERE", ".", " ", "2020."}], "  ", 
      "\[IndentingNewLine]", "Numerical", " ", "Example", " ", "2"}], ":", 
     " ", 
     RowBox[{"Agricultural", " ", "Pollution"}]}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"Quit", "[", "]"}]}]], "Input",
 CellChangeTimes->{{3.7626005245226784`*^9, 3.762600526093836*^9}, {
  3.8125441532338324`*^9, 3.8125441979767914`*^9}},
 CellLabel->
  "In[357]:=",ExpressionUUID->"851da910-5bec-44c6-af44-3b60a1f63ce3"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Functional", " ", "Forms"}], "*)"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"Unrestricted", " ", "Profit", " ", 
      RowBox[{"Function", ":", " ", 
       RowBox[{
        RowBox[{"production", " ", "f", " ", "is", " ", "cobb"}], "-", 
        RowBox[{
        "douglas", " ", "in", " ", "aggregate", " ", "nitrogen", " ", "n", 
         " ", "and", " ", "land", " ", "k"}]}], ":", " ", "f"}]}], "=", 
     RowBox[{"A", "*", 
      RowBox[{"(", 
       RowBox[{"n", "^", "\[Zeta]"}], ")"}], "*", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"k", "^", "\[Beta]"}], ")"}], ".", " ", "Note"}], " ", "that",
       " ", "parameter", " ", "symbols", " ", "here", " ", "may", " ", 
      "differ", " ", "from", " ", "the", " ", "main", " ", "text"}]}], ",", 
    " ", 
    RowBox[{
    "as", " ", "some", " ", "parameters", " ", "in", " ", "the", " ", "text", 
     " ", "were", " ", "changed", " ", "during", " ", "the", " ", "revision", 
     " ", "process"}], ",", " ", 
    RowBox[{"but", " ", "we", " ", 
     RowBox[{"didn", "'"}], "t", " ", "want", " ", "to", " ", "change", " ", 
     "them", " ", "here", " ", "for", " ", "fear", " ", "of", " ", 
     "introducing", " ", "errors", " ", "into", " ", "the", " ", 
     "program"}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"f", "=", 
     RowBox[{"A", "*", 
      RowBox[{"(", 
       RowBox[{"n", "^", "\[Zeta]"}], ")"}], "*", 
      RowBox[{"(", 
       RowBox[{"k", "^", "\[Beta]"}], ")"}]}]}], ";"}], "\[IndentingNewLine]",
    "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Unrestricted", " ", 
      RowBox[{"Profits", ".", " ", "These"}], " ", "are", " ", "denoted", " ",
       "\[Pi]x", " ", "to", " ", "indicate", " ", "that", " ", "profits", " ",
       "depend", " ", "on", " ", "input", " ", "use"}], ",", " ", 
     RowBox[{
      RowBox[{
       RowBox[{
        RowBox[{
         RowBox[{
         "which", " ", "is", " ", "typically", " ", "denoted", " ", "by", " ",
           "a", " ", "vector", " ", 
          RowBox[{"x", ".", " ", "We"}], " ", 
          RowBox[{"can", "'"}], "t", " ", "define", " ", "\[Pi]", " ", "as", 
          " ", "profits", " ", "in", " ", "Mathematica", " ", "because", " ", 
          "Mathematica", " ", "only", " ", "recognizes", " ", "\[Pi]", " ", 
          "as", " ", "3.14"}], "..."}], ".", "  ", "Output"}], " ", "price"}],
       " ", "=", " ", 
      RowBox[{
      "p", " ", "and", " ", "input", " ", "prices", " ", "are", " ", "wn", 
       " ", "and", " ", "wk"}]}]}], "*)"}], " ", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]x", "=", 
     RowBox[{
      RowBox[{"p", "*", "f"}], "-", 
      RowBox[{"wn", "*", "n"}], " ", "-", " ", 
      RowBox[{"wk", "*", "k"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
      RowBox[{"Damage", " ", 
       RowBox[{"Function", ":", " ", "Damages"}]}], " ", "=", " ", 
      RowBox[{"d1", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"total", " ", "emissions"}], ")"}], "^", "2."}], " ", 
       "Agricultural", " ", "nonpoint", " ", "source", " ", "emissions", " ", 
       "from", " ", "corn", " ", "production"}]}], ",", " ", "e", ",", " ", 
     RowBox[{
      RowBox[{
      "are", " ", "only", " ", "part", " ", "of", " ", "the", " ", "total", 
       " ", "emissions", " ", "going", " ", "into", " ", "the", " ", 
       "Chesapeake", " ", 
       RowBox[{"Bay", ".", "  ", "Denote"}], " ", "PS", " ", "emissions", " ",
        "and", " ", "other", " ", "ag", " ", "emisisons", " ", "by", " ", 
       "eOther", " ", "and", " ", "let", " ", "e"}], " ", "=", " ", 
      RowBox[{
       RowBox[{"ag", " ", "emisisons", " ", "from", " ", "corn", " ", 
        RowBox[{"production", ".", "  ", "Then"}], " ", "total", " ", 
        "emissions"}], " ", "=", " ", 
       RowBox[{"e", " ", "+", " ", 
        RowBox[{
         RowBox[{"eOther", ".", "  ", "Holding"}], " ", "eOther", " ", 
         "fixed"}]}]}]}], ",", " ", 
     RowBox[{
     "we", " ", "can", " ", "define", " ", "a", " ", "damage", " ", 
      "function", " ", "that", " ", "only", " ", "depends", " ", "on", " ", 
      "the", " ", "variable", " ", "component", " ", "e"}]}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"dam", "=", " ", 
     RowBox[{
      RowBox[{"d1", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"eOther", "+", "e"}], ")"}], "^", "2"}]}], "-", 
      RowBox[{"d1", "*", 
       RowBox[{"eOther", "^", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "First", " ", "and", " ", "second", " ", "derivatives", " ", "of", " ", 
     "the", " ", "damage", " ", "function", " ", "to", " ", "use", " ", 
     "later"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ddamde", "=", 
     RowBox[{"D", "[", 
      RowBox[{"dam", ",", "e"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"d2damde", "=", 
     RowBox[{"D", "[", 
      RowBox[{
       RowBox[{"D", "[", 
        RowBox[{"dam", ",", "e"}], "]"}], ",", "e"}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Emissions", " ", 
      RowBox[{"Function", ":", " ", "z"}]}], "=", 
     RowBox[{"\[Phi]", "*", "n"}]}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"z", "=", 
     RowBox[{"\[Phi]", "*", "n"}]}], ";"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.762601359011119*^9, 3.7626013999602137`*^9}, {
   3.7626014394821653`*^9, 3.762601460051222*^9}, {3.762603158941807*^9, 
   3.762603168027898*^9}, {3.7626032170039997`*^9, 3.7626033085568438`*^9}, {
   3.762603530754552*^9, 3.762603555344552*^9}, {3.7626041132135515`*^9, 
   3.7626044371071186`*^9}, {3.7626045621306195`*^9, 3.762604567104117*^9}, {
   3.762604602546661*^9, 3.7626046676171675`*^9}, {3.7626047030527105`*^9, 
   3.76260470774918*^9}, {3.762605020688471*^9, 3.7626050586542673`*^9}, {
   3.7626053711975183`*^9, 3.7626053729086895`*^9}, {3.7626054433377314`*^9, 
   3.762605472536651*^9}, {3.762611148100151*^9, 3.762611150603401*^9}, {
   3.7626112599353333`*^9, 3.7626112620075407`*^9}, {3.7626155364969463`*^9, 
   3.762615596297926*^9}, {3.762615704429738*^9, 3.7626157113264275`*^9}, {
   3.7626159475360465`*^9, 3.762615958823175*^9}, {3.762616058652157*^9, 
   3.7626160777400656`*^9}, 3.7626163028285723`*^9, {3.7626180494592175`*^9, 
   3.762618090836355*^9}, {3.7626184376170297`*^9, 3.762618443455613*^9}, {
   3.7626857239793053`*^9, 3.762685750576965*^9}, {3.7626859888357887`*^9, 
   3.762686009473852*^9}, {3.76268604805971*^9, 3.7626860481557198`*^9}, {
   3.7626863313020315`*^9, 3.762686341702071*^9}, {3.7626954981366234`*^9, 
   3.7626954982376337`*^9}, {3.7626973353313246`*^9, 
   3.7626973606418552`*^9}, {3.762697450931883*^9, 3.762697453948185*^9}, {
   3.762697742953082*^9, 3.762697744999287*^9}, 3.7626978407988663`*^9, 
   3.7626979206178474`*^9, {3.7627058269296203`*^9, 3.7627058435262794`*^9}, {
   3.7627065314610662`*^9, 3.76270653619954*^9}, {3.7627071110870233`*^9, 
   3.762707111283043*^9}, {3.7632205761877966`*^9, 3.76322057632181*^9}, 
   3.7632207476479406`*^9, {3.763220807035879*^9, 3.763220828137989*^9}, {
   3.7632208772919035`*^9, 3.7632208810612807`*^9}, {3.7632234406682158`*^9, 
   3.763223444358585*^9}, {3.763223493553504*^9, 3.76322350481763*^9}, {
   3.763223700258172*^9, 3.763223702945441*^9}, 3.7632238196311083`*^9, {
   3.7632259643239765`*^9, 3.763225964381976*^9}, {3.7632260988649764`*^9, 
   3.763226098973976*^9}, {3.7632261566369762`*^9, 3.7632261567389765`*^9}, {
   3.7632262107219763`*^9, 3.7632262109189763`*^9}, {3.771862015530128*^9, 
   3.77186209725495*^9}, {3.771862140845389*^9, 3.7718621867768507`*^9}, {
   3.7718623247375145`*^9, 3.7718623437487054`*^9}, {3.7718630131537437`*^9, 
   3.771863059704212*^9}, {3.771866012991952*^9, 3.771866132549163*^9}, 
   3.7718662321391773`*^9, {3.7718727370487576`*^9, 3.7718727385497723`*^9}, {
   3.77187286443104*^9, 3.771872876659163*^9}, {3.771873524222684*^9, 
   3.7718735540659847`*^9}, {3.771873788098355*^9, 3.7718737890883646`*^9}, {
   3.7719354653145747`*^9, 3.771935473690659*^9}, 3.771935647948414*^9, {
   3.7719408545768127`*^9, 3.771940855536755*^9}, {3.771940911431431*^9, 
   3.771940911802408*^9}, {3.7719470436979403`*^9, 3.771947054805052*^9}, {
   3.7719478925819407`*^9, 3.771947892771943*^9}, {3.771953408443739*^9, 
   3.7719534101947565`*^9}, 3.771953665015914*^9, {3.771955180235084*^9, 
   3.7719551803650856`*^9}, {3.772190060391234*^9, 3.7721900615682464`*^9}, {
   3.772297394122628*^9, 3.7722974278939676`*^9}, {3.772298053238286*^9, 
   3.7722980699274535`*^9}, {3.7722981211419687`*^9, 
   3.7722981291880493`*^9}, {3.7722991839836693`*^9, 
   3.7722991998328285`*^9}, {3.772302325496222*^9, 3.7723023255862217`*^9}, {
   3.7723026757978053`*^9, 3.772302701508085*^9}, {3.772302892614168*^9, 
   3.7723029038342896`*^9}, {3.805227985652342*^9, 3.805227985922124*^9}, {
   3.805280204622018*^9, 3.8052802047479477`*^9}, {3.805280258754227*^9, 
   3.8052802589890947`*^9}, {3.8052845076897583`*^9, 3.805284507838673*^9}, {
   3.8057965393409457`*^9, 3.805796641664087*^9}, {3.805796864362136*^9, 
   3.8057968889910455`*^9}, {3.8057976788809347`*^9, 3.805797861143779*^9}, {
   3.805800805467689*^9, 3.8058009157010307`*^9}, {3.805800979687319*^9, 
   3.8058009952212276`*^9}, {3.8058010277032394`*^9, 
   3.8058010391160345`*^9}, {3.805801150246728*^9, 3.8058011893392043`*^9}, {
   3.805801269761946*^9, 3.8058013079499035`*^9}, {3.805801349788344*^9, 
   3.8058013670066385`*^9}, {3.805807005102093*^9, 3.8058070277031417`*^9}, {
   3.806157617033944*^9, 3.806157630081351*^9}, {3.8061576923559113`*^9, 
   3.8061576994290953`*^9}, 3.8061577553124943`*^9, {3.8125442373096375`*^9, 
   3.8125442392272015`*^9}, {3.8125442839158816`*^9, 
   3.8125444454349165`*^9}, {3.8125444801205916`*^9, 
   3.8125445585835333`*^9}, {3.812544611390255*^9, 3.812544716672307*^9}, {
   3.8125447511272316`*^9, 3.8125448427354193`*^9}, {3.812545903757107*^9, 
   3.812545917708151*^9}, {3.8125459478880434`*^9, 3.8125459496987567`*^9}, {
   3.812545997995769*^9, 3.812546056580249*^9}, {3.8127431791514735`*^9, 
   3.8127431796067753`*^9}},
 CellLabel->"In[1]:=",ExpressionUUID->"395c95fd-1bbf-41d2-8e50-37061b2a08ce"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Compute", " ", "Restricted", " ", 
     RowBox[{"Profits", ".", " ", "The"}], " ", "value", " ", "of", " ", "n", 
     " ", "is", " ", "determined", " ", "solely", " ", "by", " ", "the", " ", 
     "emissions", " ", 
     RowBox[{"constraint", ".", " ", "Then"}], " ", "the", " ", "value", " ", 
     "of", " ", "k", " ", "is", " ", "chosen", " ", "to", " ", "maximize", 
     " ", 
     RowBox[{"profits", ".", " ", "These"}], " ", "input", " ", "responses", 
     " ", "are", " ", "defined", " ", "n1e", " ", "and", " ", "k1e"}], ",", 
    " ", 
    RowBox[{
    "with", " ", "the", " ", "resulting", " ", "restricted", " ", "profit", 
     " ", "function", " ", "defined", " ", "\[Pi]e", " ", 
     RowBox[{"(", 
      RowBox[{
      "not", " ", "to", " ", "be", " ", "confused", " ", "with", " ", "a", 
       " ", "partial", " ", "derivative", " ", "with", " ", "respect", " ", 
       "to", " ", "e"}], ")"}]}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"sol\[Pi]e", "=", 
     RowBox[{"Solve", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"e", " ", "\[Equal]", " ", "z"}], "}"}], ",", 
       RowBox[{"{", "n", "}"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"n1e", " ", "=", " ", 
     RowBox[{"Simplify", "[", 
      RowBox[{"n", "/.", 
       RowBox[{"sol\[Pi]e", "[", 
        RowBox[{"[", "1", "]"}], "]"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"k1e", "=", 
     RowBox[{"k", "/.", 
      RowBox[{
       RowBox[{"Solve", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"D", "[", 
           RowBox[{
            RowBox[{"\[Pi]x", "/.", 
             RowBox[{"n", "\[Rule]", "n1e"}]}], ",", "k"}], "]"}], "\[Equal]",
           "0"}], ",", "k"}], "]"}], "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]e", "=", 
     RowBox[{"Simplify", "[", 
      RowBox[{
       RowBox[{"\[Pi]x", "/.", 
        RowBox[{"n", "\[Rule]", "n1e"}]}], "/.", 
       RowBox[{"k", "\[Rule]", "k1e"}]}], "]"}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.8057967789615545`*^9, 3.80579681129389*^9}, {
   3.805796912536686*^9, 3.805796985117097*^9}, {3.805801326819343*^9, 
   3.8058013394130955`*^9}, {3.8058085130416136`*^9, 
   3.8058085481538568`*^9}, {3.805808580529256*^9, 3.805808609178534*^9}, {
   3.8058152889199133`*^9, 3.805815291257531*^9}, {3.8125448838500805`*^9, 
   3.812544981891383*^9}, {3.812546070003182*^9, 3.812546072347189*^9}, 
   3.812723809112219*^9, 3.8127438139806895`*^9},
 CellLabel->"In[7]:=",ExpressionUUID->"daf9159f-7516-4d32-b91a-decf05d4fe44"],

Cell[BoxData[
 TemplateBox[{
  "Solve","ifun",
   "\"Inverse functions are being used by \
\\!\\(\\*RowBox[{\\\"Solve\\\"}]\\), so some solutions may not be found; use \
Reduce for complete solution information.\"",2,9,1,28351761546097222402,
   "Local"},
  "MessageTemplate"]], "Message", "MSG",
 CellChangeTimes->{{3.8058085733923674`*^9, 3.8058086105444646`*^9}, 
   3.8058089883509617`*^9, 3.805810135961753*^9, 3.8058150544594793`*^9, 
   3.8058150893631163`*^9, 3.8058154933215046`*^9, 3.8058257117365932`*^9, 
   3.805881567941851*^9, 3.805881782249466*^9, 3.805889614038216*^9, 
   3.8058912683176785`*^9, 3.8058913629449077`*^9, 3.8058914843246326`*^9, 
   3.806064924755869*^9, 3.8061542753769493`*^9, 3.812719477128557*^9, 
   3.8127206503390684`*^9, 3.8127207862053347`*^9, 3.8127210843850365`*^9, 
   3.812722171473647*^9, 3.8127224703608246`*^9, 3.8127233251205916`*^9, 
   3.8127235351151915`*^9, 3.812723739896644*^9, 3.8127238161592283`*^9, 
   3.812724180703201*^9, 3.8127245081064324`*^9, 3.8127254456918445`*^9, 
   3.8127268478911347`*^9, {3.812743809315879*^9, 3.8127438172771144`*^9}, 
   3.8127442788316813`*^9, 3.812748190722706*^9},
 CellLabel->
  "During evaluation of \
In[7]:=",ExpressionUUID->"b8ac0353-cf74-49a5-b1ed-89d483392b61"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "Calibration", " ", "of", " ", "the", " ", "Profit", " ", "and", " ", 
    "Emissions", " ", "Function"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Parmater", " ", 
     RowBox[{"list", ":", " ", "A"}]}], ",", " ", "\[Zeta]", ",", " ", 
    "\[Beta]", ",", " ", "p", ",", " ", "wn", ",", " ", "wk", ",", " ", 
    RowBox[{"\[Phi]", " ", 
     RowBox[{"(", 
      RowBox[{"7", " ", "parameters"}], ")"}], "\[IndentingNewLine]", "The", 
     " ", "parameters", " ", "are", " ", "calibrated", " ", "using", " ", "7",
      " ", "calibration", " ", "parameters", " ", "obtained", " ", "from", 
     " ", "published", " ", 
     RowBox[{"data", ".", " ", "Some"}], " ", "of", " ", "these", " ", 
     "calibration", " ", "parameters", " ", "are", " ", "treated", " ", "as", 
     " ", "if", " ", "they", " ", "are", " ", "known", " ", "with", " ", 
     "certainty", " ", "and", " ", "the", " ", "rest", " ", "are", " ", 
     "treated", " ", "as", " ", "uncertain", " ", "but", " ", "with", " ", 
     "a", " ", "known", " ", 
     RowBox[{"distribution", ".", "\[IndentingNewLine]", "Known"}], " ", 
     "calibration", " ", 
     RowBox[{"parameters", ":", " ", "p"}]}], ",", " ", "e0", ",", " ", "k0", 
    ",", " ", 
    RowBox[{"and", " ", "TVCn"}], ",", " ", 
    RowBox[{
    "where", " ", "TVCn", " ", "is", " ", "the", " ", "total", " ", "costs", 
     " ", "associated", " ", "with", " ", "nitrogen", " ", "use", " ", 
     "\[IndentingNewLine]", "\[IndentingNewLine]", " ", "Uncertain", " ", 
     "calibration", " ", 
     RowBox[{"parameters", ":", " ", 
      RowBox[{"n0", " ", 
       RowBox[{"(", 
        RowBox[{"initial", " ", "n", " ", "use"}], ")"}]}]}]}], ",", " ", 
    RowBox[{"y0", " ", 
     RowBox[{"(", 
      RowBox[{"initial", " ", "output"}], ")"}]}], ",", " ", 
    RowBox[{"and", " ", "\[Pi]0", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{
       "initial", " ", "profits", " ", "as", " ", "a", " ", "fraction", " ", 
        "of", " ", "initial", " ", "revenue"}], ",", " ", 
       RowBox[{"i", ".", "e", "."}], ",", " ", 
       RowBox[{"\[Pi]x", "/", 
        RowBox[{"(", 
         RowBox[{"p", "*", "y0"}], ")"}]}]}], ")"}], "   ", 
     "\[IndentingNewLine]", "\[IndentingNewLine]", "Calibration", " ", 
     RowBox[{"details", ":", " ", "\[IndentingNewLine]", 
      RowBox[{"1.", " ", "p", " ", "is", " ", "determined", " ", "by", " ", 
       RowBox[{"p", ".", " ", "2."}], " ", "\[Phi]", " ", "is", " ", 
       "determined", " ", "by", " ", "e0", " ", "and", " ", 
       RowBox[{"n0", ".", " ", "3."}], " ", "wn", " ", "is", " ", 
       "determined", " ", "by", " ", "n0", " ", "and", " ", 
       RowBox[{"TVCn", ".", " ", "4."}], " ", "wk", " ", "is", " ", 
       "determined", " ", "by", " ", "\[Pi]0"}]}]}], ",", " ", "TVCn", ",", 
    " ", 
    RowBox[{
     RowBox[{"and", " ", "k0", " ", "determine", " ", 
      RowBox[{"wk", ".", " ", "6"}], " ", "and", " ", "7.", " ", "y0"}], " ", 
     "=", " ", 
     RowBox[{"f", 
      RowBox[{"(", 
       RowBox[{"n0", ",", "k0"}], ")"}], " ", "and", " ", "FOCs", " ", "with",
       " ", "respect", " ", "to", " ", "n", " ", "and", " ", "k"}]}], ",", 
    " ", 
    RowBox[{"evaluated", " ", "at", " ", "n0", " ", "and", " ", "k0"}], ",", 
    " ", 
    RowBox[{"determine", " ", "A"}], ",", " ", "\[Zeta]", ",", " ", 
    RowBox[{"and", " ", "\[Beta]"}]}], "*)"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Calibration", " ", "emissions", " ", "function"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"\[Phi]", " ", "=", 
     RowBox[{"e0", "/", "n0"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"Calibration", " ", "of", " ", "Profit", " ", "Function"}], 
    "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"wn", " ", "=", 
     RowBox[{"TVCn", "/", "n0"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"wk", " ", "=", " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"p", "*", "y0", "*", 
         RowBox[{"(", 
          RowBox[{"1", "-", "\[Pi]0"}], ")"}]}], " ", "-", "TVCn"}], ")"}], 
      "/", "k0"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Zeta]", "=", 
     RowBox[{"TVCn", "/", 
      RowBox[{"(", 
       RowBox[{"p", "*", "y0"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Beta]", "=", 
     RowBox[{"wk", "*", 
      RowBox[{"k0", "/", 
       RowBox[{"(", 
        RowBox[{"p", "*", "y0"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"A", " ", "=", " ", 
     RowBox[{"y0", "/", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"n0", "^", "\[Zeta]"}], ")"}], "*", 
        RowBox[{"(", 
         RowBox[{"k0", "^", "\[Beta]"}], ")"}]}], ")"}]}]}], 
    ";"}]}]}]], "Input",
 CellChangeTimes->{{3.7626047203654413`*^9, 3.762604782583663*^9}, {
   3.7626048306974735`*^9, 3.7626048588102846`*^9}, {3.7626050668270845`*^9, 
   3.762605085986*^9}, {3.7626053650389023`*^9, 3.7626053780132*^9}, {
   3.762605436856084*^9, 3.762605436983096*^9}, {3.762611264895829*^9, 
   3.762611301225462*^9}, {3.7626113322555647`*^9, 3.762611538106148*^9}, {
   3.762611581442481*^9, 3.762611704166752*^9}, {3.762612075132845*^9, 
   3.762612134070738*^9}, {3.7626121713524656`*^9, 3.762612514259753*^9}, {
   3.762613892711584*^9, 3.7626138963899527`*^9}, {3.762613927848098*^9, 
   3.762613931710484*^9}, {3.7626154495092487`*^9, 3.7626155256078577`*^9}, {
   3.7626156048937855`*^9, 3.7626156675010457`*^9}, {3.7626157192552204`*^9, 
   3.7626157591672115`*^9}, {3.762615990848377*^9, 3.7626160251488066`*^9}, {
   3.76261608848514*^9, 3.76261610158245*^9}, {3.7626162320894985`*^9, 
   3.7626163155648456`*^9}, {3.762616375946883*^9, 3.7626164105123396`*^9}, {
   3.7626175887611523`*^9, 3.7626175906823444`*^9}, {3.7626181037916503`*^9, 
   3.7626181057808495`*^9}, {3.762618952287492*^9, 3.7626189523865013`*^9}, {
   3.7626192624585056`*^9, 3.76261927170243*^9}, {3.762619304687728*^9, 
   3.7626193361608753`*^9}, {3.7626196454277983`*^9, 
   3.7626196459158473`*^9}, {3.7627159771865444`*^9, 
   3.7627160020090265`*^9}, {3.7627175919780073`*^9, 3.7627175926180716`*^9}, 
   3.7627177598437924`*^9, {3.763220660803257*^9, 3.763220660973274*^9}, {
   3.7632207045316296`*^9, 3.7632207046706433`*^9}, {3.7632209086540394`*^9, 
   3.7632209164858227`*^9}, {3.7632234534144907`*^9, 3.763223483529502*^9}, 
   3.7632235370338516`*^9, {3.7632235720293503`*^9, 3.7632236012612734`*^9}, {
   3.763223718145961*^9, 3.763223759691115*^9}, {3.763223845112656*^9, 
   3.763223931219266*^9}, {3.763224858067384*^9, 3.7632248586797876`*^9}, 
   3.763225680060623*^9, {3.7718626442920313`*^9, 3.7718627547341423`*^9}, {
   3.7718628263158627`*^9, 3.771862900107605*^9}, {3.7718661601794434`*^9, 
   3.7718661954598026`*^9}, {3.7718662405122623`*^9, 
   3.7718662553414135`*^9}, {3.771866362951518*^9, 3.7718664014409175`*^9}, {
   3.7718665116820593`*^9, 3.7718665431333847`*^9}, {3.7718715453515205`*^9, 
   3.7718716053141236`*^9}, {3.771871642890503*^9, 3.7718716430005035`*^9}, {
   3.7718717043241205`*^9, 3.771871741413494*^9}, {3.7718718403294897`*^9, 
   3.771871924488336*^9}, {3.7718719545516386`*^9, 3.771871981989915*^9}, {
   3.7718723601057186`*^9, 3.771872406540186*^9}, {3.771953697941267*^9, 
   3.771953750723831*^9}, {3.771954891376176*^9, 3.77195489672123*^9}, {
   3.7719551116273937`*^9, 3.7719551774740562`*^9}, {3.7721254593087983`*^9, 
   3.772125461430402*^9}, {3.772190108157728*^9, 3.7721902128578157`*^9}, {
   3.772190787888874*^9, 3.772190790548901*^9}, {3.772191261975896*^9, 
   3.7721912720660024`*^9}, {3.77219130474736*^9, 3.7721913092274075`*^9}, {
   3.7721913450267987`*^9, 3.7721913617899785`*^9}, {3.7722980663254175`*^9, 
   3.772298067507429*^9}, {3.772298114497902*^9, 3.772298148691246*^9}, {
   3.7723618820613537`*^9, 3.7723618937024755`*^9}, {3.8049813849151616`*^9, 
   3.8049813926078744`*^9}, {3.8052215028728385`*^9, 
   3.8052215132880588`*^9}, {3.805221560964285*^9, 3.8052215634645844`*^9}, {
   3.8052844125155706`*^9, 3.805284426565645*^9}, {3.8052844654706497`*^9, 
   3.805284496981041*^9}, 3.8053239159280777`*^9, {3.8053240105523396`*^9, 
   3.8053240119072943`*^9}, {3.8053240815434093`*^9, 
   3.8053240852642903`*^9}, {3.8053257313988075`*^9, 3.805325740023769*^9}, 
   3.8053262709702187`*^9, {3.8054566402352977`*^9, 3.8054566403082542`*^9}, {
   3.8054584689358077`*^9, 3.805458469366558*^9}, {3.8057209581202135`*^9, 
   3.8057209664514217`*^9}, 3.8057429207991753`*^9, {3.805797663360433*^9, 
   3.805797667757662*^9}, {3.805797919652852*^9, 3.8057979458552885`*^9}, {
   3.8057979780247645`*^9, 3.805798133776072*^9}, {3.8057981663735294`*^9, 
   3.805798267040193*^9}, {3.8057983423055487`*^9, 3.805798404032685*^9}, {
   3.805799147726745*^9, 3.8057991776389894`*^9}, {3.8057992133994265`*^9, 
   3.8057992304140587`*^9}, {3.8057992857812214`*^9, 3.805799417664786*^9}, {
   3.805799453037318*^9, 3.805799551340742*^9}, {3.805799620304222*^9, 
   3.8057996393391085`*^9}, {3.8057997583114014`*^9, 
   3.8057997816163845`*^9}, {3.805799835201795*^9, 3.8057999821381693`*^9}, {
   3.8058000127505016`*^9, 3.805800018162442*^9}, {3.8058000497903757`*^9, 
   3.8058001653861723`*^9}, {3.805800196918562*^9, 3.8058001996265755`*^9}, {
   3.80580058150287*^9, 3.8058005840983396`*^9}, {3.8058010464839487`*^9, 
   3.8058010842923937`*^9}, {3.805801222901251*^9, 3.8058012585047836`*^9}, {
   3.8058013576681385`*^9, 3.805801478675071*^9}, {3.8058015716046762`*^9, 
   3.8058016658152113`*^9}, {3.805801717934466*^9, 3.8058018464768105`*^9}, {
   3.805801885745576*^9, 3.805801889692298*^9}, {3.805801956525302*^9, 
   3.8058019985234947`*^9}, 3.805802126097386*^9, {3.8058070482565217`*^9, 
   3.805807105690911*^9}, {3.8058071712440476`*^9, 3.8058071930020733`*^9}, {
   3.8058083293651276`*^9, 3.8058084352082167`*^9}, {3.8058086256231976`*^9, 
   3.8058086280968575`*^9}, 3.80580869030869*^9, {3.80580906199099*^9, 
   3.805809065741844*^9}, {3.8058233738887234`*^9, 3.8058233740516295`*^9}, {
   3.8058255873336477`*^9, 3.805825635703712*^9}, {3.8061576371461225`*^9, 
   3.806157658228054*^9}, {3.806157705543284*^9, 3.8061577251934247`*^9}, {
   3.8125453093904324`*^9, 3.8125453211224566`*^9}, {3.8125453562992887`*^9, 
   3.8125455571204824`*^9}, {3.812545654705533*^9, 3.8125457006468287`*^9}, {
   3.8125461364122295`*^9, 3.8125461728752375`*^9}, {3.8125467615518403`*^9, 
   3.812546797518343*^9}, {3.8125468370384226`*^9, 3.8125468418898335`*^9}, {
   3.812547496406344*^9, 3.8125475081265616`*^9}, {3.8125476703891497`*^9, 
   3.8125476711264453`*^9}, {3.8125477093430495`*^9, 
   3.8125477103101864`*^9}, {3.8125478730587907`*^9, 3.812547873948574*^9}, {
   3.8127241681730366`*^9, 3.8127241707796884`*^9}, {3.812724217591119*^9, 
   3.812724219859177*^9}},
 CellLabel->"In[11]:=",ExpressionUUID->"7b96dfab-41a9-4b2b-a3f4-5f1c44b9c1db"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Paramater", " ", "Values", " ", "and", " ", 
     RowBox[{"Distributions", ".", " ", "The"}], " ", "distributions", " ", 
     "have", " ", "been", " ", "chosen", " ", "to", " ", "ensure", " ", 
     "combinations", " ", "of", " ", "realized", " ", "uncertain", " ", 
     "calibration", " ", "parameters", " ", "do", " ", "not", " ", "result", 
     " ", "in", " ", "unrealistic", " ", "calibrations"}], ",", " ", 
    RowBox[{"e", ".", "g", "."}], ",", " ", 
    RowBox[{"negative", " ", "profits"}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Known", " ", "calibration", " ", "parameters"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"e0", " ", "=", " ", 
     RowBox[{"60.6", "*", 
      RowBox[{"10", "^", "6"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"k0", "=", "1450000"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"p", " ", "=", " ", "4"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"TVCn", " ", "=", " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"4", "/", "6"}], ")"}], "*", "p", "*", "y0"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "Uncertain", " ", "calibration", " ", "parameters", " ", "are", " ", 
     "assumed", " ", "uniformly", " ", 
     RowBox[{"distributed", ".", "  ", "Here"}], " ", "we", " ", "set", " ", 
     "the", " ", "upper", " ", "and", " ", "lower", " ", "bounds", " ", "for",
      " ", "these", " ", "distributions"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]0high", "=", "0.35"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]0low", "=", "0.03"}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"n0high", "=", 
     RowBox[{"220", "*", "k0"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"n0low", "=", 
     RowBox[{"150", "*", "k0"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"y0low", "=", 
     RowBox[{"120", "*", "k0"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"y0high", "=", 
     RowBox[{"198", "*", "k0"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Specify", " ", "Damage", " ", "Function", " ", "parameters"}], 
     ",", " ", 
     RowBox[{"d1", " ", "and", " ", "eOther"}], ",", " ", 
     RowBox[{"which", " ", "are", " ", "assumed", " ", "known"}]}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"e0", "=", 
     RowBox[{"60.6", "*", 
      RowBox[{"10", "^", "6"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"eOther", "=", 
     RowBox[{"5.45", "*", 
      RowBox[{"10", "^", "6"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"d1", "=", 
     RowBox[{"4.14", "*", 
      RowBox[{"10", "^", 
       RowBox[{"-", "8"}]}]}]}], ";"}], "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.771945645436839*^9, 3.7719456759851465`*^9}, {
   3.771945765842052*^9, 3.771945780223197*^9}, {3.7719458438538365`*^9, 
   3.7719459529969344`*^9}, {3.771945987121278*^9, 3.7719460231696424`*^9}, {
   3.771946143650855*^9, 3.771946157992999*^9}, {3.7719462056504793`*^9, 
   3.7719462196286197`*^9}, {3.7719464857902975`*^9, 
   3.7719464871713114`*^9}, {3.7719465286847305`*^9, 3.771946541275857*^9}, {
   3.771946589206339*^9, 3.7719465905773525`*^9}, {3.7719466768972216`*^9, 
   3.7719466825002775`*^9}, {3.771947101017517*^9, 3.771947116017668*^9}, {
   3.771947462554614*^9, 3.771947463224621*^9}, {3.771953774094083*^9, 
   3.7719537831351805`*^9}, {3.77195654583317*^9, 3.7719565938046527`*^9}, {
   3.7719567646853714`*^9, 3.771956922001958*^9}, {3.7719569954206963`*^9, 
   3.7719570445931907`*^9}, {3.771957184887602*^9, 3.771957211985874*^9}, {
   3.771957466020444*^9, 3.771957481099595*^9}, {3.771957737296174*^9, 
   3.7719580308441267`*^9}, {3.7719580828076496`*^9, 3.77195815042233*^9}, {
   3.7720999365849175`*^9, 3.7720999479573374`*^9}, {3.772100029483081*^9, 
   3.772100111492425*^9}, {3.7721001529572973`*^9, 3.7721004399822016`*^9}, {
   3.7721004773286676`*^9, 3.7721005612256145`*^9}, {3.772100765897974*^9, 
   3.772100803946441*^9}, 3.7721008704961576`*^9, {3.7721022443125706`*^9, 
   3.772102292750656*^9}, {3.772102342389943*^9, 3.7721023461495495`*^9}, {
   3.772103468212721*^9, 3.7721035080239906`*^9}, {3.7721035808449183`*^9, 
   3.7721035955089445`*^9}, {3.772103941330352*^9, 3.7721039606899853`*^9}, {
   3.772104323674423*^9, 3.7721043275588303`*^9}, {3.7721043630674925`*^9, 
   3.772104394595148*^9}, {3.772104618988942*^9, 3.772104619207343*^9}, {
   3.772105316703168*^9, 3.7721053187779713`*^9}, {3.7721288862856784`*^9, 
   3.7721288893432837`*^9}, {3.772128939297572*^9, 3.7721289473003855`*^9}, {
   3.772190256252266*^9, 3.772190269032401*^9}, {3.7721910274314237`*^9, 
   3.772191072471899*^9}, {3.7721911068942623`*^9, 3.772191108194276*^9}, {
   3.7721911567547884`*^9, 3.772191158464805*^9}, {3.77219209945481*^9, 
   3.772192109034911*^9}, {3.7721935547424607`*^9, 3.772193557112489*^9}, {
   3.772193670113718*^9, 3.77219367938482*^9}, {3.7721942685691004`*^9, 
   3.7721943044634614`*^9}, {3.7721944039804635`*^9, 3.772194486226299*^9}, {
   3.7722014482391434`*^9, 3.7722014491791525`*^9}, {3.7722960550761433`*^9, 
   3.772296056317156*^9}, {3.772296118779784*^9, 3.7722961202517986`*^9}, {
   3.772296219686799*^9, 3.772296221998822*^9}, {3.772296591580544*^9, 
   3.7722966566161995`*^9}, {3.772297408631774*^9, 3.7722974204088926`*^9}, 
   3.7723013247959237`*^9, {3.772302832592514*^9, 3.772302882793061*^9}, {
   3.7723030641260347`*^9, 3.772303065206045*^9}, {3.7723033382999516`*^9, 
   3.77230335504012*^9}, {3.7723041740827503`*^9, 3.772304178687797*^9}, {
   3.7723042779097977`*^9, 3.7723042781698*^9}, {3.772361790312385*^9, 
   3.772361810611599*^9}, {3.7723618464029756`*^9, 3.7723618485629997`*^9}, {
   3.7723619281028385`*^9, 3.772361928672844*^9}, {3.772364508114134*^9, 
   3.772364513704194*^9}, {3.7723650471090364`*^9, 3.772365063900218*^9}, {
   3.772366347928517*^9, 3.7723663501373672`*^9}, {3.772366473969943*^9, 
   3.772366479541692*^9}, {3.772369241388362*^9, 3.772369272016683*^9}, {
   3.772371365622365*^9, 3.772371367554384*^9}, {3.772371857721321*^9, 
   3.772371861774362*^9}, {3.7723723279080524`*^9, 3.77237232967807*^9}, {
   3.782145107335912*^9, 3.78214510851007*^9}, {3.782145479004854*^9, 
   3.7821454984665895`*^9}, {3.782145538938282*^9, 3.78214554579613*^9}, {
   3.7821455759244523`*^9, 3.7821456365699625`*^9}, {3.804982918530241*^9, 
   3.8049829217625704`*^9}, {3.8049832205796776`*^9, 3.80498322218786*^9}, 
   3.8052196088341894`*^9, {3.805219851631589*^9, 3.8052198884257765`*^9}, {
   3.8052200494570312`*^9, 3.805220057031231*^9}, {3.8052210377959657`*^9, 
   3.8052210413314333`*^9}, {3.805221619227539*^9, 3.8052216267200427`*^9}, {
   3.805280020062905*^9, 3.805280022652321*^9}, {3.8052802764093094`*^9, 
   3.8052802782560034`*^9}, {3.805284585096345*^9, 3.805284587532608*^9}, {
   3.8052846551680126`*^9, 3.805284657068558*^9}, {3.8052880802078133`*^9, 
   3.8052880808295774`*^9}, {3.80532310479282*^9, 3.805323105284681*^9}, {
   3.8053231742745066`*^9, 3.8053232018269367`*^9}, {3.805323281434894*^9, 
   3.805323318661711*^9}, {3.805323625973091*^9, 3.805323692355037*^9}, {
   3.8053238479429955`*^9, 3.8053238576565275`*^9}, {3.8053246405915833`*^9, 
   3.8053246412250786`*^9}, {3.8053247355420823`*^9, 
   3.8053247989626255`*^9}, {3.805324969355548*^9, 3.805324971130247*^9}, {
   3.80532505791402*^9, 3.805325059130088*^9}, 3.8053254969260497`*^9, {
   3.80532564642485*^9, 3.805325647536193*^9}, {3.805325852154505*^9, 
   3.8053258545561347`*^9}, {3.8053259331428533`*^9, 
   3.8053259670271206`*^9}, {3.805326063758665*^9, 3.805326081761572*^9}, {
   3.8053261990147357`*^9, 3.805326200290702*^9}, {3.8054557078918805`*^9, 
   3.8054557176962137`*^9}, 3.8057430241115828`*^9, {3.8057432340200877`*^9, 
   3.8057432423299913`*^9}, {3.805743764510635*^9, 3.805743779477542*^9}, {
   3.805809152039339*^9, 3.805809153098726*^9}, {3.805809280177262*^9, 
   3.805809294380577*^9}, 3.8058101119238653`*^9, {3.805810636890899*^9, 
   3.8058106460815945`*^9}, {3.8058151064973507`*^9, 
   3.8058151071310225`*^9}, {3.8058154369427958`*^9, 
   3.8058154511522245`*^9}, {3.8058235096415405`*^9, 
   3.8058235364500837`*^9}, {3.8058255335589485`*^9, 
   3.8058255513206005`*^9}, {3.8058258425321355`*^9, 
   3.8058258433777385`*^9}, {3.8058259243068385`*^9, 
   3.8058259262323103`*^9}, {3.8058260498812904`*^9, 3.80582605164123*^9}, {
   3.805826113336261*^9, 3.8058261155860386`*^9}, {3.8058261893326654`*^9, 
   3.805826191872037*^9}, {3.805826312367351*^9, 3.805826337035124*^9}, {
   3.805826494650182*^9, 3.8058264955066957`*^9}, {3.805826581713108*^9, 
   3.805826581755065*^9}, {3.805826660137314*^9, 3.8058266805396757`*^9}, {
   3.805826745159343*^9, 3.8058267452023067`*^9}, {3.805881681293738*^9, 
   3.8058816820993114`*^9}, {3.805881740245278*^9, 3.805881749384099*^9}, {
   3.8058818194745655`*^9, 3.80588182466741*^9}, 3.80588191826019*^9, {
   3.8058819740720973`*^9, 3.805881975437641*^9}, {3.805891249091342*^9, 
   3.8058912497913413`*^9}, {3.8058913394239397`*^9, 3.8058913461300063`*^9}, 
   3.805891431135112*^9, {3.805891463157134*^9, 3.805891464799027*^9}, 
   3.806157730372625*^9, {3.8125471629127884`*^9, 3.812547192438036*^9}, {
   3.8125472352030945`*^9, 3.8125473519619565`*^9}, 3.8125473880005074`*^9, {
   3.8125474412132225`*^9, 3.8125474867835045`*^9}, {3.8125475674541492`*^9, 
   3.8125475705267987`*^9}, {3.81254764694098*^9, 3.8125476595960846`*^9}, {
   3.8125477263725066`*^9, 3.8125478685455074`*^9}, {3.8127194859939675`*^9, 
   3.812719487265628*^9}, {3.8127210688520346`*^9, 3.81272107009264*^9}, {
   3.8127254248086443`*^9, 3.812725426102031*^9}, {3.81272547477307*^9, 
   3.81272547490176*^9}, {3.8127255201474366`*^9, 3.812725522173805*^9}, {
   3.8127257655213013`*^9, 3.8127257658016043`*^9}, {3.8127268585694885`*^9, 
   3.8127268944078712`*^9}, {3.812726928527226*^9, 3.812726988987383*^9}, {
   3.812727072702164*^9, 3.812727203525241*^9}, {3.812744253970953*^9, 
   3.8127442618459835`*^9}, {3.8127481744958735`*^9, 3.812748177831024*^9}, 
   3.8127482493258657`*^9},
 CellLabel->"In[32]:=",ExpressionUUID->"71d7e32a-4ae8-4ebb-b902-48a504476939"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "Calculate", " ", "expected", " ", "restricted", " ", "profit", " ", 
     RowBox[{"function", ".", " ", "Expected"}], " ", "values", " ", "will", 
     " ", "be", " ", "calculated", " ", "by", " ", "integrating", " ", "over",
      " ", "the", " ", "independent"}], ",", " ", 
    RowBox[{"uncertain", " ", "parameters", " ", "n0"}], ",", " ", "\[Pi]0", 
    ",", " ", 
    RowBox[{"and", " ", 
     RowBox[{"y0", ".", " ", "Integrations"}], " ", "will", " ", "be", " ", 
     "performed", " ", "using", " ", 
     RowBox[{"Simpson", "'"}], "s", " ", "rule"}], ",", " ", 
    RowBox[{
     RowBox[{"where", " ", "'"}], 
     RowBox[{"v", "'"}], " ", "denotes", " ", "the", " ", "number", " ", "of",
      " ", "intervals", " ", "used", " ", "for", " ", "this", " ", 
     RowBox[{"rule", ".", " ", "We"}], " ", "use", " ", 
     RowBox[{"Simpson", "'"}], "s", " ", "rule", " ", "so", " ", "that", " ", 
     "we", " ", "have", " ", "an", " ", "analytical", " ", "expected", " ", 
     "profit", " ", "function", " ", "that", " ", "we", " ", "can", " ", 
     RowBox[{"differentiate", ".", " ", "An"}], " ", "alternative", " ", 
     "approach", " ", "might", " ", "have", " ", "been", " ", "to", " ", 
     "use", " ", "the", " ", "integration", " ", "command", " ", 
     RowBox[{"Integrate", "[", "]"}]}], ",", " ", 
    RowBox[{
    "which", " ", "we", " ", "used", " ", "for", " ", "Example", " ", "1"}], 
    ",", " ", 
    RowBox[{
    "but", " ", "the", " ", "model", " ", "is", " ", "too", " ", "nonlinear", 
     " ", "and", " ", "has", " ", "too", " ", "many", " ", "uncertain", " ", 
     "variables", " ", "for", " ", "that", " ", "to", " ", "work"}]}], "*)"}],
   "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
   "Specify", " ", "weights", " ", "for", " ", "the", " ", "uncertain", " ", 
    "variables"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"weight\[Pi]", "=", 
     RowBox[{"1", "/", 
      RowBox[{"(", 
       RowBox[{"\[Pi]0high", "-", "\[Pi]0low"}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"weighty", " ", "=", " ", 
     RowBox[{"1", "/", 
      RowBox[{"(", 
       RowBox[{"y0high", " ", "-", " ", "y0low"}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"weightn", "=", 
     RowBox[{"1", "/", 
      RowBox[{"(", 
       RowBox[{"n0high", "-", "n0low"}], ")"}]}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"Number", " ", "of", " ", "intervals", " ", "for", " ", 
     RowBox[{"Simpson", "'"}], "s", " ", "rule"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"v", " ", "=", "20"}], ";"}], "  ", "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Simpson", "'"}], "s", " ", "Rule", " ", "Procedure", " ", "to", 
     " ", "Calculate", " ", "expected", " ", "restricted", " ", "profits"}], 
    "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{
     RowBox[{"SimpsonsSum", "[", 
      RowBox[{"g_", ",", "num_", ",", "a_", ",", "b_"}], "]"}], ":=", 
     RowBox[{"Module", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"width", "=", 
          RowBox[{
           RowBox[{"(", 
            RowBox[{"b", "-", "a"}], ")"}], "/", "num"}]}], ",", 
         "coefficient"}], "}"}], ",", 
       RowBox[{
        RowBox[{
         RowBox[{"coefficient", "[", 
          RowBox[{"i_", "?", "EvenQ"}], "]"}], "=", "2"}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"coefficient", "[", 
          RowBox[{"i_", "?", "OddQ"}], "]"}], "=", "4"}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"N", "[", 
         RowBox[{
          RowBox[{"width", "/", "3"}], "*", 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"g", "[", "a", "]"}], "+", 
            RowBox[{"g", "[", "b", "]"}], "+", 
            RowBox[{"Sum", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"coefficient", "[", "i", "]"}], "*", 
               RowBox[{"g", "[", 
                RowBox[{"a", "+", 
                 RowBox[{"i", "*", "width"}]}], "]"}]}], ",", 
              RowBox[{"{", 
               RowBox[{"i", ",", "1", ",", 
                RowBox[{"num", "-", "1"}]}], "}"}]}], "]"}]}], ")"}]}], 
         "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]ey", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"\[Pi]e", "/.", 
        RowBox[{"y0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]ey", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Pi]ey", ",", "v", ",", "y0low", ",", "y0high"}], "]"}], "*",
       "weighty"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Variables", "[", "E\[Pi]ey", "]"}], ";"}], "\[IndentingNewLine]",
    "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]e\[Pi]", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Pi]ey", "/.", 
        RowBox[{"\[Pi]0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]e\[Pi]", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Pi]e\[Pi]", ",", "v", ",", "\[Pi]0low", ",", "\[Pi]0high"}],
        "]"}], "*", "weight\[Pi]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Variables", "[", "E\[Pi]e\[Pi]", "]"}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]en", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Pi]e\[Pi]", "/.", 
        RowBox[{"n0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]e", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Pi]en", ",", "v", ",", "n0low", ",", "n0high"}], "]"}], "*",
       "weightn"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Variables", "[", "E\[Pi]e", "]"}], ";"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.7626164261148996`*^9, 3.762616431479436*^9}, {
   3.762616960763359*^9, 3.7626169846517477`*^9}, 3.762620613725619*^9, {
   3.7626217084280777`*^9, 3.7626219571579485`*^9}, {3.7626219991911516`*^9, 
   3.7626220029435263`*^9}, {3.762622928537077*^9, 3.7626229299762206`*^9}, {
   3.7626229631535378`*^9, 3.762622973105533*^9}, {3.762623194074628*^9, 
   3.7626232082260427`*^9}, {3.7626233765068693`*^9, 
   3.7626234179300113`*^9}, {3.7626236092341394`*^9, 
   3.7626237771509295`*^9}, {3.7626238688681*^9, 3.762623873114525*^9}, 
   3.7626240736675787`*^9, {3.7626241947186823`*^9, 3.7626242618083906`*^9}, {
   3.7626252350447044`*^9, 3.762625240412241*^9}, {3.762626111556347*^9, 
   3.762626115321723*^9}, {3.762626165884779*^9, 3.762626188373028*^9}, {
   3.7626263222354126`*^9, 3.7626263287860675`*^9}, {3.762626509506138*^9, 
   3.762626512760463*^9}, {3.762626660728258*^9, 3.7626267282160063`*^9}, {
   3.7626267692501097`*^9, 3.7626268039075747`*^9}, {3.762626982564439*^9, 
   3.762626991066289*^9}, {3.7626270253797197`*^9, 3.762627102766458*^9}, {
   3.762627298717051*^9, 3.762627299245104*^9}, {3.7626954355533657`*^9, 
   3.7626954384596562`*^9}, {3.7626955557143803`*^9, 3.762695557326542*^9}, {
   3.762797704213747*^9, 3.7627977404453697`*^9}, {3.763216149348157*^9, 
   3.7632161494591675`*^9}, {3.7632161837105923`*^9, 3.763216183812603*^9}, {
   3.7722974538812294`*^9, 3.772297496170655*^9}, {3.7722975398800945`*^9, 
   3.7722975446431427`*^9}, {3.7722979625373735`*^9, 
   3.7722980307130594`*^9}, {3.7722983915966897`*^9, 
   3.7722986585923753`*^9}, {3.772300525371566*^9, 3.772300632651724*^9}, {
   3.772300675975108*^9, 3.7723006760750556`*^9}, {3.772300727995536*^9, 
   3.7723007896128764`*^9}, {3.772301815812448*^9, 3.7723018709612217`*^9}, {
   3.7723019129202223`*^9, 3.7723019675312223`*^9}, 3.772302014952222*^9, {
   3.7723020809722223`*^9, 3.772302157750222*^9}, {3.772302392504222*^9, 
   3.772302395306222*^9}, {3.772361707410512*^9, 3.772361760411069*^9}, {
   3.772361946643033*^9, 3.772362012103726*^9}, {3.81254794073449*^9, 
   3.812548164889605*^9}, {3.812548200337165*^9, 3.8125483519812365`*^9}, {
   3.812584962240368*^9, 3.8125850864551105`*^9}, {3.812722271761804*^9, 
   3.812722272365241*^9}, 3.8127438364120026`*^9},
 CellLabel->"In[31]:=",ExpressionUUID->"cbf0001c-ea66-4e0d-ac41-d76bb60a4987"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
   "Compute", " ", "the", " ", "Optimal", " ", "Standard", " ", "using", " ", 
    "the", " ", "Nonlinear", " ", "Model"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"E\[Pi]e0", "=", 
     RowBox[{"E\[Pi]e", "/.", 
      RowBox[{"e", "\[Rule]", "e0"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Edam0", "=", 
     RowBox[{"dam", "/.", 
      RowBox[{"e", "\[Rule]", "e0"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ESNB", "=", " ", 
     RowBox[{"(", 
      RowBox[{"E\[Pi]e", "-", "dam"}], ")"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"solehat", "=", 
     RowBox[{"FindRoot", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{"ESNB", ",", "e"}], "]"}], "\[Equal]", "0"}], ",", 
       RowBox[{"{", 
        RowBox[{"e", ",", 
         RowBox[{"0.5", "*", "e0"}]}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ehat", "=", 
     RowBox[{"e", "/.", 
      RowBox[{"solehat", "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "Economic", " ", "values", " ", "associated", " ", "with", " ", "the", 
     " ", "optimal", " ", "standard"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]ehat", "=", 
     RowBox[{"E\[Pi]e", "/.", 
      RowBox[{"e", "\[Rule]", "ehat"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Edamhat", "=", 
     RowBox[{"dam", "/.", 
      RowBox[{"e", "\[Rule]", "ehat"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ESNBhat", "=", 
     RowBox[{"ESNB", "/.", 
      RowBox[{"e", "\[Rule]", "ehat"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "Marginal", " ", "damages", " ", "at", " ", "the", " ", "optimal", " ", 
     "standard"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"M\[Pi]ehat", "=", 
     RowBox[{
      RowBox[{"D", "[", 
       RowBox[{"E\[Pi]e", ",", "e"}], "]"}], "/.", 
      RowBox[{"e", "\[Rule]", "ehat"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "Sign", " ", "of", " ", "Relative", " ", "Slopes", " ", "at", " ", "the", 
     " ", "Optimal", " ", "Standard"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"relativeslopes", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", 
           RowBox[{"D", "[", 
            RowBox[{
             RowBox[{"D", "[", 
              RowBox[{"E\[Pi]e", ",", "e"}], "]"}], ",", "e"}], "]"}]}], "-", 
          RowBox[{"D", "[", 
           RowBox[{
            RowBox[{"D", "[", 
             RowBox[{"dam", ",", "e"}], "]"}], ",", "e"}], "]"}]}], ")"}], 
        "/", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"D", "[", 
           RowBox[{"dam", ",", "e"}], "]"}], ",", "e"}], "]"}]}], ")"}], "/.", 
      RowBox[{"e", "\[Rule]", "ehat"}]}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{
  3.762621856042838*^9, {3.762624086778889*^9, 3.762624135266738*^9}, {
   3.762625294485648*^9, 3.7626253171729164`*^9}, {3.7626253517363725`*^9, 
   3.762625356391838*^9}, {3.762625412025401*^9, 3.7626254121124096`*^9}, {
   3.7626261949906893`*^9, 3.762626296698859*^9}, {3.7626271629704776`*^9, 
   3.7626271699361744`*^9}, {3.7626274166118393`*^9, 
   3.7626274384420223`*^9}, {3.7626902868445463`*^9, 3.762690297835645*^9}, 
   3.772302307834222*^9, {3.772302432374222*^9, 3.772302471115574*^9}, {
   3.7723025058169527`*^9, 3.772302616077154*^9}, {3.7723029830441513`*^9, 
   3.7723029834041553`*^9}, {3.7723031282857327`*^9, 
   3.7723031284357343`*^9}, {3.7723036045156307`*^9, 
   3.7723037960035563`*^9}, {3.7723044205872355`*^9, 3.772304463289666*^9}, 
   3.772307837643253*^9, {3.7723079326786194`*^9, 3.7723079343166227`*^9}, {
   3.7723079662342787`*^9, 3.772307993066326*^9}, {3.772308093124902*^9, 
   3.772308144168191*^9}, {3.7723081777238503`*^9, 3.772308215928317*^9}, {
   3.772308813632767*^9, 3.7723088255355883`*^9}, {3.772362441418254*^9, 
   3.7723624929688005`*^9}, {3.805222946660144*^9, 3.805222970119108*^9}, {
   3.805229334791565*^9, 3.8052293491024857`*^9}, {3.805229675826316*^9, 
   3.8052296771701756`*^9}, {3.8053233934784603`*^9, 
   3.8053234184961987`*^9}, {3.8053239692214403`*^9, 3.805323989805935*^9}, {
   3.8053241609688625`*^9, 3.805324174530509*^9}, {3.805815170833398*^9, 
   3.8058152191646967`*^9}, {3.8058817963387513`*^9, 3.805881801808539*^9}, {
   3.805890025644539*^9, 3.805890034655382*^9}, {3.8125483733725986`*^9, 
   3.812548375246929*^9}, {3.8125838161101623`*^9, 3.8125839160597143`*^9}, {
   3.812584719876359*^9, 3.8125848040965433`*^9}, {3.812584868825792*^9, 
   3.812584911696044*^9}, {3.812720729495115*^9, 3.812720730014053*^9}, {
   3.812720990358919*^9, 3.812720991935655*^9}, {3.812721136089105*^9, 
   3.8127211373253098`*^9}},
 CellLabel->"In[45]:=",ExpressionUUID->"1a7c9eee-15a8-4798-a586-1912b00536eb"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "Develop", " ", "quadratic", " ", "approimxations", " ", "around", " ", 
     RowBox[{"ehat", ".", " ", "This"}], " ", "follows", " ", 
     RowBox[{"Weitzman", "'"}], "s", " ", "approach"}], ",", " ", 
    RowBox[{
    "except", " ", "that", " ", "we", " ", "allow", " ", "the", " ", "slope", 
     " ", "of", " ", "marginal", " ", "profits", " ", "to", " ", "be", " ", 
     RowBox[{"uncertain", ".", " ", "Use"}], " ", "this", " ", "model", " ", 
     "to", " ", "Compute", " ", "the", " ", "TRUE", " ", "Optimal", " ", 
     "Tax", " ", "and", " ", "associated", " ", "TRUE", " ", "measure", " ", 
     "of", " ", "ESNB"}]}], "*)"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "Approximate", " ", "profits", " ", "around", " ", "the", " ", "optimal", 
     " ", "standard", " ", 
     RowBox[{"ehat", ".", "  ", "Here"}], " ", "we", " ", "need", " ", "to", 
     " ", "approximate", " ", "the", " ", "quadratic", " ", "coefficients"}], 
    ",", " ", "\[Pi]0", ",", " ", "\[Pi]e1", ",", " ", 
    RowBox[{"and", " ", "\[Pi]e2", " ", 
     RowBox[{"(", 
      RowBox[{
      "with", " ", "\[Pi]e1", " ", "and", " ", "\[Pi]e2", " ", "being", " ", 
       "the", " ", "first", " ", "and", " ", "second", " ", "derivatives", 
       " ", "of", " ", "restricted", " ", "profits", " ", "with", " ", 
       "respect", " ", "to", " ", "e"}], ")"}]}], ",", " ", 
    RowBox[{
    "each", " ", "of", " ", "which", " ", "we", " ", "treat", " ", "as", " ", 
     RowBox[{"uncertain", ".", " ", "Note"}], " ", "that"}], ",", " ", 
    RowBox[{"for", " ", "now"}], ",", " ", 
    RowBox[{"we", " ", 
     RowBox[{"don", "'"}], "t", " ", "bother", " ", "trying", " ", "to", " ", 
     "decompose", " ", "\[Pi]e1", " ", "and", " ", "\[Pi]e2", " ", "as", " ", 
     "the", " ", "combination", " ", "of", " ", "an", " ", "expected", " ", 
     "value", " ", "and", " ", "a", " ", "random", " ", "term", " ", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"e", ".", "g", "."}], ",", " ", 
        RowBox[{"\[Pi]e1", " ", "=", " ", 
         RowBox[{"E\[Pi]e1", "*", 
          RowBox[{"(", 
           RowBox[{"1", "+", "\[Alpha]"}], ")"}], " ", "in", " ", "the", " ", 
          "text"}]}], ",", " ", 
        RowBox[{
        "but", " ", "here", " ", "we", " ", "simply", " ", "define", " ", 
         "the", " ", "term", " ", "\[Pi]e1", " ", "without", " ", "computing",
          " ", "E\[Pi]e1", " ", "and", " ", "\[Alpha]", " ", "separately"}]}],
        ")"}], ".", " ", "There"}], " ", "is", " ", "no", " ", "need", " ", 
     "for", " ", "such", " ", "a", " ", "decomposition", " ", "in", " ", 
     "order", " ", "to", " ", "produce", " ", "numerical", " ", "results"}]}],
    "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"\[Pi]e0", "=", 
     RowBox[{"\[Pi]e", "/.", 
      RowBox[{"e", "\[Rule]", "ehat"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]e1", "=", 
     RowBox[{
      RowBox[{"D", "[", 
       RowBox[{"\[Pi]e", ",", "e"}], "]"}], "/.", 
      RowBox[{"e", "\[Rule]", "ehat"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]e2", "=", 
     RowBox[{
      RowBox[{"D", "[", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{"\[Pi]e", ",", "e"}], "]"}], ",", "e"}], "]"}], "/.", 
      RowBox[{"e", "\[Rule]", "ehat"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]equad", " ", "=", " ", 
     RowBox[{"\[Pi]e0", " ", "+", 
      RowBox[{"\[Pi]e1", "*", 
       RowBox[{"(", 
        RowBox[{"e", "-", "ehat"}], ")"}]}], "+", 
      RowBox[{"0.5", "*", "\[Pi]e2", "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"e", "-", "ehat"}], ")"}], "^", "2"}]}]}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Compute", " ", "the", " ", 
      RowBox[{"firm", "'"}], "s", " ", "response", " ", "to", " ", "the", " ",
       "tax"}], ",", " ", 
     RowBox[{"i", ".", "e", "."}], ",", " ", 
     RowBox[{"its", " ", "response", " ", "function"}], ",", " ", 
     RowBox[{
      RowBox[{
       RowBox[{"based", " ", "on", " ", "the", " ", "FOC", " ", "d", 
        RowBox[{
         RowBox[{"(", "\[Pi]equad", ")"}], "/", "de"}]}], " ", "-", " ", 
       "\[Tau]"}], " ", "=", " ", "0"}]}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"e\[Tau]", "=", 
     RowBox[{"Simplify", "[", 
      RowBox[{"ehat", "+", 
       RowBox[{"\[Tau]", "/", "\[Pi]e2"}], "-", 
       RowBox[{"\[Pi]e1", "/", "\[Pi]e2"}]}], "]"}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.81272124028568*^9, 3.812721275581077*^9}, {
  3.8127439258610096`*^9, 3.8127439263880944`*^9}, {3.812744111895283*^9, 
  3.8127441563440137`*^9}, {3.8127467769075565`*^9, 
  3.8127468337006607`*^9}},ExpressionUUID->"23f8ba2f-2397-4804-9220-\
95de22826350"],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Define", " ", "the", " ", "optimal"}], ",", " ", 
    RowBox[{
    "approximated", " ", "profit", " ", "and", " ", "damage", " ", 
     "functions"}], ",", " ", 
    RowBox[{
    "and", " ", "associated", " ", "social", " ", "net", " ", "benefits"}], 
    ",", " ", 
    RowBox[{"given", " ", "the", " ", 
     RowBox[{"firm", "'"}], "s", " ", "response", " ", "function"}]}], "*)"}],
   "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"\[Pi]equad\[Tau]", "=", 
     RowBox[{"\[Pi]equad", "/.", 
      RowBox[{"e", "\[Rule]", "e\[Tau]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"damquad\[Tau]", "=", 
     RowBox[{"dam", "/.", 
      RowBox[{"e", "\[Rule]", "e\[Tau]"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"snb\[Tau]", "=", 
     RowBox[{"\[Pi]equad\[Tau]", "-", "damquad\[Tau]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Variables", "[", "\[Pi]equad\[Tau]", "]"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Variables", "[", "damquad\[Tau]", "]"}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"Use", " ", 
     RowBox[{"Simpson", "'"}], "s", " ", "Rule", " ", "to", " ", "calculate", 
     " ", "ESNB", " ", "for", " ", "a", " ", "given", " ", "\[Tau]"}], "*)"}],
    "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"snb\[Tau]y", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"snb\[Tau]", "/.", 
        RowBox[{"y0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Esnb\[Tau]y", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"snb\[Tau]y", ",", "v", ",", "y0low", ",", "y0high"}], "]"}], 
      "*", "weighty"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"E\[Pi]e\[Theta]1", "=", 
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Pi]1\[Theta]", ",", "v", ",", "a1", ",", "b1"}], "]"}]}], 
     ";"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Variables", "[", "Esnb\[Tau]y", "]"}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"snb\[Tau]\[Pi]", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Esnb\[Tau]y", "/.", 
        RowBox[{"\[Pi]0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Esnb\[Tau]\[Pi]", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{
       "snb\[Tau]\[Pi]", ",", "v", ",", "\[Pi]0low", ",", "\[Pi]0high"}], 
       "]"}], "*", "weight\[Pi]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Variables", "[", "Esnb\[Tau]\[Pi]", "]"}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"snb\[Tau]n", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Esnb\[Tau]\[Pi]", "/.", 
        RowBox[{"n0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Esnb\[Tau]", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"snb\[Tau]n", ",", "v", ",", "n0low", ",", "n0high"}], "]"}], 
      "*", "weightn"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Variables", "[", "Esnb\[Tau]", "]"}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "Solve", " ", "for", " ", "the", " ", "Truly", " ", "Optimal", " ", 
     "Tax"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Sol\[Tau]", "=", 
     RowBox[{"FindRoot", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"D", "[", 
          RowBox[{"Esnb\[Tau]", ",", "\[Tau]"}], "]"}], "\[Equal]", "0"}], 
        "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Tau]", ",", "4.2"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Tau]1", "=", 
     RowBox[{"\[Tau]", "/.", 
      RowBox[{"Sol\[Tau]", "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "Compute", " ", "the", " ", "Truly", " ", "optimal", " ", "ESNB", " ", 
     "and", " ", "true", " ", "\[CapitalDelta]"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Esnb\[Tau]", "/.", 
     RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[CapitalDelta]true", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Esnb\[Tau]", "/.", 
        RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1"}]}], ")"}], "-", "ESNBhat"}]}],
     ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Compute", " ", "the", " ", "risk", " ", 
      RowBox[{"term", " ", "'"}], 
      RowBox[{"r", "'"}]}], ":", " ", 
     RowBox[{
      RowBox[{"if", " ", "r"}], " ", ">", " ", 
      RowBox[{
      "0", " ", "then", " ", "the", " ", "intercept", " ", "and", " ", 
       "slope", " ", "are", " ", "negatively", " ", "correlated"}]}]}], 
    "*)"}], " ", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"risk", "=", 
     RowBox[{"\[Tau]1", "-", "M\[Pi]ehat"}]}], ";"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.772304954217104*^9, 3.7723049579701414`*^9}, {
   3.772305000729572*^9, 3.772305072665295*^9}, {3.7723616064694443`*^9, 
   3.772361683799262*^9}, {3.7723617207106504`*^9, 3.77236172814073*^9}, {
   3.7723620809244504`*^9, 3.7723621188848543`*^9}, {3.772362227796003*^9, 
   3.772362260226343*^9}, {3.7723623522073154`*^9, 3.7723623817176247`*^9}, {
   3.772371208986788*^9, 3.7723712480341825`*^9}, {3.7723726761185617`*^9, 
   3.7723727272130795`*^9}, {3.772372796959781*^9, 3.7723727970707827`*^9}, {
   3.7821461649217477`*^9, 3.7821461829326477`*^9}, {3.812585132025096*^9, 
   3.812585229904625*^9}, {3.8125853142341986`*^9, 3.8125856872477407`*^9}, {
   3.812585799562511*^9, 3.8125859119272976`*^9}, {3.8125859542318015`*^9, 
   3.8125859635021048`*^9}, {3.812586016501753*^9, 3.8125860445814915`*^9}, {
   3.8125861133013296`*^9, 3.8125861668361373`*^9}, {3.8125862002011623`*^9, 
   3.812586269281098*^9}, 3.8127212185299473`*^9},
 CellLabel->"In[60]:=",ExpressionUUID->"8f646553-56c0-4a66-9846-7ab75b07e623"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Construct", " ", 
     RowBox[{"Weitzman", "'"}], "s", " ", "analysis", " ", "based", " ", "on",
      " ", "a", " ", "quadratic", " ", "approximation", " ", "of", " ", "the",
      " ", "profit", " ", "function"}], ",", " ", 
    RowBox[{
     RowBox[{
     "given", " ", "ehat", " ", "and", " ", "assuming", " ", "\[Gamma]"}], 
     "=", 
     RowBox[{
     "1.", " ", "This", " ", "latter", " ", "assumption", " ", "simply", " ", 
      "means", " ", "substituting", " ", "E", 
      RowBox[{"{", "\[Pi]e2", "}"}], " ", "in", " ", "for", " ", "\[Pi]e2", 
      " ", "in", " ", "the", " ", "quadratic", " ", "profit", " ", 
      "function"}]}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"\[Pi]equadW", " ", "=", " ", 
     RowBox[{"\[Pi]e0", " ", "+", 
      RowBox[{"\[Pi]e1", "*", 
       RowBox[{"(", 
        RowBox[{"e", "-", "ehat"}], ")"}]}], "+", 
      RowBox[{"0.5", "*", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"D", "[", 
          RowBox[{
           RowBox[{"D", "[", 
            RowBox[{"E\[Pi]e", ",", "e"}], "]"}], ",", "e"}], "]"}], "/.", 
         RowBox[{"e", "\[Rule]", " ", "ehat"}]}], ")"}], "*", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"e", "-", "ehat"}], ")"}], "^", "2"}]}]}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"The", " ", 
     RowBox[{"firm", "'"}], "s", " ", "predicted", " ", "response", " ", "to",
      " ", "the", " ", "tax", " ", "under", " ", 
     RowBox[{"Weitzman", "'"}], "s", " ", "assumptions"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"sole\[Tau]W", "=", 
     RowBox[{"Solve", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{"\[Pi]equadW", ",", "e"}], "]"}], "\[Equal]", "\[Tau]"}], 
       ",", "e"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"e\[Tau]W", "=", 
     RowBox[{"Simplify", "[", 
      RowBox[{"e", "/.", 
       RowBox[{"sole\[Tau]W", "[", 
        RowBox[{"[", "1", "]"}], "]"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{"Expected", " ", "SNB", " ", "under", " ", "Weitzman"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]equad\[Tau]W", "=", 
     RowBox[{"\[Pi]equadW", "/.", 
      RowBox[{"e", "\[Rule]", "e\[Tau]W"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"damquad\[Tau]W", "=", 
     RowBox[{"dam", "/.", 
      RowBox[{"e", "\[Rule]", "e\[Tau]W"}]}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"snb\[Tau]W", "=", 
     RowBox[{"\[Pi]equad\[Tau]W", "-", "damquad\[Tau]W"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Variables", "[", "\[Pi]equad\[Tau]W", "]"}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Variables", "[", "damquad\[Tau]W", "]"}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"snb\[Tau]yW", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"snb\[Tau]W", "/.", 
        RowBox[{"y0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Esnb\[Tau]yW", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"snb\[Tau]yW", ",", "v", ",", "y0low", ",", "y0high"}], "]"}], 
      "*", "weighty"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"E\[Pi]e\[Theta]1", "=", 
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Pi]1\[Theta]", ",", "v", ",", "a1", ",", "b1"}], "]"}]}], 
     ";"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Variables", "[", "Esnb\[Tau]yW", "]"}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"snb\[Tau]\[Pi]W", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Esnb\[Tau]yW", "/.", 
        RowBox[{"\[Pi]0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Esnb\[Tau]\[Pi]W", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{
       "snb\[Tau]\[Pi]W", ",", "v", ",", "\[Pi]0low", ",", "\[Pi]0high"}], 
       "]"}], "*", "weight\[Pi]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Variables", "[", "Esnb\[Tau]\[Pi]W", "]"}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"snb\[Tau]nW", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Esnb\[Tau]\[Pi]W", "/.", 
        RowBox[{"n0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Esnb\[Tau]W", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"snb\[Tau]nW", ",", "v", ",", "n0low", ",", "n0high"}], "]"}], 
      "*", "weightn"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Variables", "[", "Esnb\[Tau]W", "]"}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
      RowBox[{"Weitzman", "'"}], "s", " ", "Optimal", " ", "Tax", " ", 
      "Calculation"}], ";", " ", 
     RowBox[{
      RowBox[{"OPTIMUM", " ", "SHOULD", " ", "BE", " ", "TAX"}], " ", "=", 
      " ", 
      RowBox[{"MD", " ", "FROM", " ", "OPTIMAL", " ", "STANDARD"}]}]}], 
    "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Sol\[Tau]W", "=", 
     RowBox[{"FindRoot", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"D", "[", 
          RowBox[{"Esnb\[Tau]W", ",", "\[Tau]"}], "]"}], "\[Equal]", "0"}], 
        "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Tau]", ",", "M\[Pi]ehat"}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Tau]1W", "=", 
     RowBox[{"\[Tau]", "/.", 
      RowBox[{"Sol\[Tau]W", "[", 
       RowBox[{"[", "1", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
    "Expected", " ", "Social", " ", "Net", " ", "Benefits", " ", "and", " ", 
     "\[CapitalDelta]"}], "*)"}], 
   RowBox[{
    RowBox[{"Esnb\[Tau]W", "/.", 
     RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1W"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[CapitalDelta]W", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"Esnb\[Tau]W", "/.", 
        RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1W"}]}], ")"}], "-", 
      "ESNBhat"}]}], ";"}], "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.7723121872376385`*^9, 3.77231225667336*^9}, {
  3.772312326109082*^9, 3.772312491110572*^9}, {3.7723593840121555`*^9, 
  3.772359391932237*^9}, {3.772359462082954*^9, 3.7723595033743763`*^9}, {
  3.7723596384157677`*^9, 3.772359647354861*^9}, {3.7723599953924627`*^9, 
  3.772360006676578*^9}, {3.772360597707811*^9, 3.772360664828517*^9}, {
  3.7723607368022766`*^9, 3.7723607965409064`*^9}, {3.7723615677580376`*^9, 
  3.772361583330202*^9}, {3.7723620912975597`*^9, 3.77236213365501*^9}, {
  3.772362197445683*^9, 3.772362222125943*^9}, {3.812586296550744*^9, 
  3.8125863675705376`*^9}, {3.812588699931773*^9, 3.812588740411727*^9}},
 CellLabel->"In[79]:=",ExpressionUUID->"ea5b2a10-3df4-4037-83d4-434ad2e8e09a"],

Cell[CellGroupData[{

Cell[BoxData[{"\"\<Summary of Results\>\"", "\[IndentingNewLine]", \
"\"\<Unregulated ESNB\>\"", "\[IndentingNewLine]", 
 RowBox[{"ESNB0", "=", 
  RowBox[{"Esnb\[Tau]", "/.", 
   RowBox[{
   "\[Tau]", "\[Rule]", 
    "0"}]}]}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", "\"\
\<Optimal Standard as a proportion of initial emissions\>\"", "\
\[IndentingNewLine]", 
 RowBox[{"ehat", "/", 
  "e0"}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", \
"\"\<ESNB with Standard\>\"", "\[IndentingNewLine]", 
 RowBox[{"ESNBhat", "=", 
  RowBox[{"ESNB", "/.", 
   RowBox[{
   "e", "\[Rule]", 
    "ehat"}]}]}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", \
"\"\<Expected Gain in SNB from regulating under the standard\>\"", "\
\[IndentingNewLine]", 
 RowBox[{"ESNBhat", "-", 
  "ESNB0"}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", \
"\"\<Marginal damages\>\"", "\[IndentingNewLine]", 
 RowBox[{"M\[Pi]ehat", "=", 
  RowBox[{
   RowBox[{"D", "[", 
    RowBox[{"E\[Pi]e", ",", "e"}], "]"}], "/.", 
   RowBox[{
   "e", "\[Rule]", 
    "ehat"}]}]}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", \
"\"\<Sign of Relative Slopes\>\"", "\[IndentingNewLine]", 
 RowBox[{"relativeslopes", "=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", 
        RowBox[{"D", "[", 
         RowBox[{
          RowBox[{"D", "[", 
           RowBox[{"E\[Pi]e", ",", "e"}], "]"}], ",", "e"}], "]"}]}], "-", 
       RowBox[{"D", "[", 
        RowBox[{
         RowBox[{"D", "[", 
          RowBox[{"dam", ",", "e"}], "]"}], ",", "e"}], "]"}]}], ")"}], "/", 
     RowBox[{
      RowBox[{"D", "[", 
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{"E\[Pi]e", ",", "e"}], "]"}], ",", "e"}], "]"}], "^", 
      "2"}]}], ")"}], "/.", 
   RowBox[{
   "e", "\[Rule]", 
    "e0"}]}]}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", \
"\"\<Slope of MD at optimal standard\>\"", "\[IndentingNewLine]", 
 RowBox[{"d2damde", "/.", 
  RowBox[{
  "e", "\[Rule]", " ", 
   "ehat"}]}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", "\"\
\<Truly Optimal Tax\>\"", "\[IndentingNewLine]", 
 RowBox[{"\[Tau]1", "=", 
  RowBox[{"\[Tau]", "/.", 
   RowBox[{"Sol\[Tau]", "[", 
    RowBox[{"[", "1", "]"}], 
    "]"}]}]}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", 
 RowBox[{"\"\<risk term 'r': if positive then intercept and slope are \
negatively correlated\>\"", " "}], "\[IndentingNewLine]", 
 RowBox[{"risk", "=", 
  RowBox[{
  "\[Tau]1", "-", 
   "M\[Pi]ehat"}]}], "\[IndentingNewLine]", "\"\< \>\"", \
"\[IndentingNewLine]", "\"\<Truly optimal ESNB\>\"", "\[IndentingNewLine]", 
 RowBox[{"Esnb\[Tau]", "/.", 
  RowBox[{
  "\[Tau]", "\[Rule]", 
   "\[Tau]1"}]}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", \
"\"\<Expected Gain in SNB from regulating under the true tax\>\"", "\
\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(", 
   RowBox[{"Esnb\[Tau]", "/.", 
    RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1"}]}], ")"}], "-", 
  "ESNB0"}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", \
"\"\<True \[CapitalDelta]\>\"", "\[IndentingNewLine]", 
 RowBox[{"\[CapitalDelta]true", "=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"Esnb\[Tau]", "/.", 
     RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1"}]}], ")"}], "-", 
   "ESNBhat"}]}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", \
"\"\< \>\"", "\[IndentingNewLine]", "\"\<Expected SNB for an unregulated \
system, assuming \[Gamma]=1\>\"", "\[IndentingNewLine]", 
 RowBox[{"ESNBW0", "=", 
  RowBox[{"(", 
   RowBox[{"Esnb\[Tau]W", "/.", 
    RowBox[{"\[Tau]", "\[Rule]", "0"}]}], 
   ")"}]}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", \
"\"\<Weitzman's Optimal Tax\>\"", "\[IndentingNewLine]", 
 RowBox[{"\[Tau]1W", "=", 
  RowBox[{"\[Tau]", "/.", 
   RowBox[{"Sol\[Tau]W", "[", 
    RowBox[{"[", "1", "]"}], 
    "]"}]}]}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", "\"\
\<W's ESNB\>\"", "\[IndentingNewLine]", 
 RowBox[{"Esnb\[Tau]W", "/.", 
  RowBox[{
  "\[Tau]", "\[Rule]", 
   "\[Tau]1W"}]}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", \
"\"\<Expected Gain in SNB from regulating under Weitzman's tax and assuming \
\[Gamma]=1 when deriving unregulated ESNB\>\"", "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(", 
   RowBox[{"Esnb\[Tau]W", "/.", 
    RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1W"}]}], ")"}], "-", 
  "ESNBW0"}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", \
"\"\<W's \[CapitalDelta]\>\"", "\[IndentingNewLine]", 
 RowBox[{"\[CapitalDelta]W", "=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"Esnb\[Tau]W", "/.", 
     RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1W"}]}], ")"}], "-", 
   "ESNBhat"}]}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", \
"\"\< \>\"", "\[IndentingNewLine]", "\"\<Correct welfare prediction from \
Weitzman's tax\>\"", "\[IndentingNewLine]", 
 RowBox[{"Esnb\[Tau]", "/.", 
  RowBox[{
  "\[Tau]", "\[Rule]", 
   "\[Tau]1W"}]}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", \
"\"\<Correct \[CapitalDelta]\>\"", "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(", 
   RowBox[{"Esnb\[Tau]", "/.", 
    RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1W"}]}], ")"}], "-", 
  "ESNBhat"}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", "\"\
\<True Expected Gain in SNB from regulating under Weitzman's tax\>\"", "\
\[IndentingNewLine]", 
 RowBox[{"Egaintrue", "=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"Esnb\[Tau]", "/.", 
     RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1W"}]}], ")"}], "-", 
   "ESNB0"}]}], "\[IndentingNewLine]", "\"\< \>\"", "\[IndentingNewLine]", \
"\"\<% diff in ESNB from policy based on correlated uncertainty\>\"", "\
\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"100", "*", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"Max", "[", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"Esnb\[Tau]", "/.", 
          RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1"}]}], ")"}], ",", "ESNBhat"}],
        "]"}], "-", 
      RowBox[{"(", 
       RowBox[{"Esnb\[Tau]", "/.", 
        RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1W"}]}], ")"}]}], ")"}], "/", 
    RowBox[{"Abs", "[", 
     RowBox[{"Max", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"Esnb\[Tau]", "/.", 
         RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1"}]}], ")"}], ",", "ESNBhat"}], 
      "]"}], "]"}]}]}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", "\"\<% diff in gain in ESNB \
from policy based on correlated uncertainty\>\"", "\[IndentingNewLine]", 
 RowBox[{"100", "*", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{
     RowBox[{"Max", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"Esnb\[Tau]", "/.", 
         RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1"}]}], ")"}], ",", "ESNBhat"}], 
      "]"}], "-", 
     RowBox[{"(", 
      RowBox[{"Esnb\[Tau]", "/.", 
       RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1W"}]}], ")"}]}], ")"}], "/", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"Abs", "[", 
      RowBox[{"Max", "[", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"Esnb\[Tau]", "/.", 
          RowBox[{"\[Tau]", "\[Rule]", "\[Tau]1"}]}], ")"}], ",", "ESNBhat"}],
        "]"}], "]"}], "-", "ESNB0"}], ")"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]"}]}]}], "Input",
 CellChangeTimes->{{3.772362521519099*^9, 3.7723625837797585`*^9}, {
   3.7723626634816055`*^9, 3.7723627141611443`*^9}, 3.772363077365005*^9, {
   3.7723723485302596`*^9, 3.7723723765785418`*^9}, {3.7723725751325455`*^9, 
   3.772372583768632*^9}, {3.772382025902175*^9, 3.77238203432726*^9}, 
   3.77238207608668*^9, {3.812584644807001*^9, 3.8125846518965387`*^9}, {
   3.812589318309933*^9, 3.81258931915977*^9}, {3.8125893550495777`*^9, 
   3.8125893589945*^9}, {3.812717795558884*^9, 3.8127177965704517`*^9}, {
   3.8127178863962717`*^9, 3.81271788838841*^9}, {3.812717924638464*^9, 
   3.8127179574229217`*^9}, {3.8127180162528315`*^9, 3.812718345755763*^9}, {
   3.8127183875954647`*^9, 3.812718429919154*^9}, {3.81274455357098*^9, 
   3.812744553922767*^9}, {3.8127446305784807`*^9, 3.812744633338146*^9}, {
   3.812744727909055*^9, 3.8127447307932844`*^9}},
 CellLabel->
  "In[234]:=",ExpressionUUID->"9a82675f-d59c-43ec-99e2-4b1b09762af4"],

Cell[BoxData["\<\"Summary of Results\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732315246*^9},
 CellLabel->
  "Out[234]=",ExpressionUUID->"7423f4b1-0b2c-4169-82f0-4ec57d895ff3"],

Cell[BoxData["\<\"Unregulated ESNB\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.81274473232123*^9},
 CellLabel->
  "Out[235]=",ExpressionUUID->"0c6a1c50-313e-419b-a82c-58af9a48579b"],

Cell[BoxData["6.628253838899659`*^7"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447323232245`*^9},
 CellLabel->
  "Out[236]=",ExpressionUUID->"16b0601f-0310-4a4e-975d-440d4537ffc9"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447323242216`*^9},
 CellLabel->
  "Out[237]=",ExpressionUUID->"d2b04fa4-20af-46c4-bfe0-005335eee82f"],

Cell[BoxData["\<\"Optimal Standard as a proportion of initial \
emissions\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447323262167`*^9},
 CellLabel->
  "Out[238]=",ExpressionUUID->"cd19888a-8180-4c46-aa0b-b6cde8cfb859"],

Cell[BoxData["0.3801222802218071`"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732330206*^9},
 CellLabel->
  "Out[239]=",ExpressionUUID->"42092a27-8fb1-4b10-8c63-e3ed52c27d87"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447323312035`*^9},
 CellLabel->
  "Out[240]=",ExpressionUUID->"4c44ea21-4065-4b39-b0fd-07f64174236d"],

Cell[BoxData["\<\"ESNB with Standard\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447323331976`*^9},
 CellLabel->
  "Out[241]=",ExpressionUUID->"c0f13d81-c851-4863-9857-8f2957d2c7c1"],

Cell[BoxData["1.073373588250863`*^8"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732366109*^9},
 CellLabel->
  "Out[242]=",ExpressionUUID->"9f4dc2a3-9d73-45f8-a786-bf7f2c728f2c"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732368104*^9},
 CellLabel->
  "Out[243]=",ExpressionUUID->"2df40f5c-4a1d-4dc0-9d58-d88a62e987c8"],

Cell[BoxData["\<\"Expected Gain in SNB from regulating under the \
standard\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447323710957`*^9},
 CellLabel->
  "Out[244]=",ExpressionUUID->"ce566f49-dc23-40fc-9ccf-7d417c0f2a92"],

Cell[BoxData["4.10548204360897`*^7"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447323730907`*^9},
 CellLabel->
  "Out[245]=",ExpressionUUID->"2c91a6a1-ba51-490f-abb9-433eaa74e5d9"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732375085*^9},
 CellLabel->
  "Out[246]=",ExpressionUUID->"10be7a1f-c009-4b1e-8753-a96232f0a8bd"],

Cell[BoxData["\<\"Marginal damages\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.81274473237708*^9},
 CellLabel->
  "Out[247]=",ExpressionUUID->"b674b698-2da9-4dba-89b0-a34c9b946fcf"],

Cell[BoxData["2.358591963023357`"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447323790746`*^9},
 CellLabel->
  "Out[248]=",ExpressionUUID->"6a5a2ffe-bd7b-47db-a82f-b5b7b5a2ccf1"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732381069*^9},
 CellLabel->
  "Out[249]=",ExpressionUUID->"d6061b84-bfab-4d2e-972a-1c85c5a33594"],

Cell[BoxData["\<\"Sign of Relative Slopes\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732383064*^9},
 CellLabel->
  "Out[250]=",ExpressionUUID->"84c0e03c-22cb-43f3-84ff-c16d735063d1"],

Cell[BoxData[
 RowBox[{"-", "3.729164241503559`*^7"}]], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732458016*^9},
 CellLabel->
  "Out[251]=",ExpressionUUID->"3db2f32c-5c86-4bcc-a592-538263016639"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732460038*^9},
 CellLabel->
  "Out[252]=",ExpressionUUID->"98b3ecb8-a03e-4c41-96d7-1fe52fdfa3e1"],

Cell[BoxData["\<\"Slope of MD at optimal standard\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447324620056`*^9},
 CellLabel->
  "Out[253]=",ExpressionUUID->"293c2cb0-fcb0-4dce-b369-406e4ea70790"],

Cell[BoxData["8.28`*^-8"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447324639997`*^9},
 CellLabel->
  "Out[254]=",ExpressionUUID->"22544025-0796-45e2-a9a4-45b7bf9ce576"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447324660225`*^9},
 CellLabel->
  "Out[255]=",ExpressionUUID->"33ed285a-78bb-4ce2-9826-fac27fdac0ef"],

Cell[BoxData["\<\"Truly Optimal Tax\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447324680166`*^9},
 CellLabel->
  "Out[256]=",ExpressionUUID->"1aeaa1b7-e00d-47db-807b-01fb29d11da0"],

Cell[BoxData["1.5486160742573527`"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447324689865`*^9},
 CellLabel->
  "Out[257]=",ExpressionUUID->"39d8c8e2-1bc4-4d44-b902-ac9bbe1ed20b"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732470981*^9},
 CellLabel->
  "Out[258]=",ExpressionUUID->"4a3ac2f6-c934-4bf2-b198-f9c4f6acede0"],

Cell[BoxData["\<\"risk term 'r': if positive then intercept and slope are \
negatively correlated\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447324729753`*^9},
 CellLabel->
  "Out[259]=",ExpressionUUID->"5fc91c3f-8abd-477d-ba78-4192d35e7eb2"],

Cell[BoxData[
 RowBox[{"-", "0.8099758887660045`"}]], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.81274473247497*^9},
 CellLabel->
  "Out[260]=",ExpressionUUID->"31d895ef-e65d-47bf-a50d-53dcb76d6036"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447324769645`*^9},
 CellLabel->
  "Out[261]=",ExpressionUUID->"2a4ead87-757e-4e7e-9ca0-c893d558f881"],

Cell[BoxData["\<\"Truly optimal ESNB\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447324801893`*^9},
 CellLabel->
  "Out[262]=",ExpressionUUID->"598fbaf1-4ec0-4008-ac8f-c1394206183e"],

Cell[BoxData["1.0463641114351752`*^8"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447324821835`*^9},
 CellLabel->
  "Out[263]=",ExpressionUUID->"bf5aea58-8446-4715-89c5-a0d1e5989cbd"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447324841785`*^9},
 CellLabel->
  "Out[264]=",ExpressionUUID->"d4e1a51b-8b30-40fb-b473-a98c80a0e592"],

Cell[BoxData["\<\"Expected Gain in SNB from regulating under the true \
tax\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447324871707`*^9},
 CellLabel->
  "Out[265]=",ExpressionUUID->"c9ffc715-cc93-446d-8c6f-45a67a00d250"],

Cell[BoxData["3.835387275452093`*^7"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447324891653`*^9},
 CellLabel->
  "Out[266]=",ExpressionUUID->"a81caffd-af12-4f8d-a8cc-61acfb99eaef"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.81274473249116*^9},
 CellLabel->
  "Out[267]=",ExpressionUUID->"0978c7b0-870c-44da-94e6-354033661617"],

Cell[BoxData["\<\"True \[CapitalDelta]\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447324931545`*^9},
 CellLabel->
  "Out[268]=",ExpressionUUID->"b7684e05-7437-4101-89b0-013867be5357"],

Cell[BoxData[
 RowBox[{"-", "2.7009476815687716`*^6"}]], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732496146*^9},
 CellLabel->
  "Out[269]=",ExpressionUUID->"88fe9d5d-c381-459e-8f50-977b402a0920"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732498141*^9},
 CellLabel->
  "Out[270]=",ExpressionUUID->"e13a8f71-755b-47e2-9a6a-5cbd3c90dcec"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325001354`*^9},
 CellLabel->
  "Out[271]=",ExpressionUUID->"bddb10ef-6eba-4142-bdd1-161741289484"],

Cell[BoxData["\<\"Expected SNB for an unregulated system, assuming \[Gamma]=1\
\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325021305`*^9},
 CellLabel->
  "Out[272]=",ExpressionUUID->"34a2ee8b-d8b5-4483-9ba2-7b37830629fd"],

Cell[BoxData["6.933137690450406`*^7"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732504125*^9},
 CellLabel->
  "Out[273]=",ExpressionUUID->"13a7ab4e-e840-4c2e-a6db-44f650e1add2"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325061197`*^9},
 CellLabel->
  "Out[274]=",ExpressionUUID->"22679df2-3f69-4121-aad1-16fa555915d1"],

Cell[BoxData["\<\"Weitzman's Optimal Tax\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732508114*^9},
 CellLabel->
  "Out[275]=",ExpressionUUID->"995b8406-4a7b-4f4d-b0a1-1a14e364696e"],

Cell[BoxData["2.358591963023357`"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.81274473251213*^9},
 CellLabel->
  "Out[276]=",ExpressionUUID->"4514676c-be1f-43dd-ab07-0b8fdf8cb745"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325141277`*^9},
 CellLabel->
  "Out[277]=",ExpressionUUID->"b5557ed8-8252-4ff9-8f82-52f3cd06ed5e"],

Cell[BoxData["\<\"W's ESNB\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325150957`*^9},
 CellLabel->
  "Out[278]=",ExpressionUUID->"811e8c7e-2750-4d39-b0a0-9030ef15523f"],

Cell[BoxData["1.088658457835091`*^8"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325171185`*^9},
 CellLabel->
  "Out[279]=",ExpressionUUID->"c6a5603d-75e4-4322-84be-395c6d1ebf07"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732519086*^9},
 CellLabel->
  "Out[280]=",ExpressionUUID->"189ab2bc-102e-4d56-b4ac-0541b23883cf"],

Cell[BoxData["\<\"Expected Gain in SNB from regulating under Weitzman's tax \
and assuming \[Gamma]=1 when deriving unregulated ESNB\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325211077`*^9},
 CellLabel->
  "Out[281]=",ExpressionUUID->"2f9759df-b04a-4072-8cb4-6c83fb4eb671"],

Cell[BoxData["3.9534468879005045`*^7"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325231028`*^9},
 CellLabel->
  "Out[282]=",ExpressionUUID->"fe0bd304-f75e-4089-9d18-9ae8266cf8da"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732525107*^9},
 CellLabel->
  "Out[283]=",ExpressionUUID->"1751448f-97f4-4bbe-9ce6-7ff470dd96ff"],

Cell[BoxData["\<\"W's \[CapitalDelta]\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325270643`*^9},
 CellLabel->
  "Out[284]=",ExpressionUUID->"a6a27a59-05a5-4d94-a8b7-6c5602a41c41"],

Cell[BoxData["1.5284869584228098`*^6"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732529087*^9},
 CellLabel->
  "Out[285]=",ExpressionUUID->"8ac3ab23-c685-49dc-ba11-48859476cc02"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325310817`*^9},
 CellLabel->
  "Out[286]=",ExpressionUUID->"7f40b1dd-b8d2-4a48-becb-25c77a04e709"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325330815`*^9},
 CellLabel->
  "Out[287]=",ExpressionUUID->"244f61c9-f2a9-4cf6-835d-0b3d63ed0f72"],

Cell[BoxData["\<\"Correct welfare prediction from Weitzman's tax\"\>"], \
"Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325350714`*^9},
 CellLabel->
  "Out[288]=",ExpressionUUID->"ac00df29-b43d-4af8-9791-7905764d6173"],

Cell[BoxData["9.414422759284021`*^7"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732537057*^9},
 CellLabel->
  "Out[289]=",ExpressionUUID->"30354b1b-8066-4d60-9e61-75232f5736cb"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732539033*^9},
 CellLabel->
  "Out[290]=",ExpressionUUID->"06aefbb0-befe-4091-b03e-4afbd4579a93"],

Cell[BoxData["\<\"Correct \[CapitalDelta]\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732541026*^9},
 CellLabel->
  "Out[291]=",ExpressionUUID->"a0c3ca86-63c9-4ef5-8159-e209f3d1f76c"],

Cell[BoxData[
 RowBox[{"-", "1.3193131232246086`*^7"}]], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325430613`*^9},
 CellLabel->
  "Out[292]=",ExpressionUUID->"2bd99fcb-b5a5-4aab-b064-7818f76cc0d4"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325460424`*^9},
 CellLabel->
  "Out[293]=",ExpressionUUID->"dec9bfbf-399e-4e94-b63d-10f41c8dfad2"],

Cell[BoxData["\<\"True Expected Gain in SNB from regulating under Weitzman's \
tax\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325480404`*^9},
 CellLabel->
  "Out[294]=",ExpressionUUID->"e120990f-d8ff-46b7-bc5c-7720d87511f6"],

Cell[BoxData["2.7861689203843616`*^7"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325493803`*^9},
 CellLabel->
  "Out[295]=",ExpressionUUID->"8f84a9cf-4759-4a27-a3c2-c44650725801"],

Cell[BoxData["\<\" \"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325513887`*^9},
 CellLabel->
  "Out[296]=",ExpressionUUID->"8582eecf-4b95-4f34-b44e-38f3c865f194"],

Cell[BoxData["\<\"% diff in ESNB from policy based on correlated \
uncertainty\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325533533`*^9},
 CellLabel->
  "Out[297]=",ExpressionUUID->"d563a08f-1f0e-4455-b259-652a071f8b2a"],

Cell[BoxData["12.291276193729725`"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732555349*^9},
 CellLabel->
  "Out[298]=",ExpressionUUID->"82d910d7-27ed-4562-8fda-ce6634e2447e"],

Cell[BoxData["\<\"% diff in gain in ESNB from policy based on correlated \
uncertainty\"\>"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.812744732559338*^9},
 CellLabel->
  "Out[299]=",ExpressionUUID->"7535d366-2349-41ae-8e6e-00c9b4b807a3"],

Cell[BoxData["32.135401134646095`"], "Output",
 CellChangeTimes->{
  3.772362849193582*^9, 3.772363079085022*^9, 3.772365117361798*^9, 
   3.7723663948730936`*^9, 3.7723693930849056`*^9, 3.772371412255839*^9, 
   3.77237191806693*^9, {3.7723723604573793`*^9, 3.772372377521551*^9}, 
   3.772372478247568*^9, 3.7723725843886385`*^9, 3.772382035038267*^9, 
   3.772382079248712*^9, 3.7723829534015183`*^9, 3.782144269738433*^9, 
   3.782145393498144*^9, 3.783255552392755*^9, 3.805743091563511*^9, 
   3.8058958450464563`*^9, 3.806064953291202*^9, 3.8061543096268625`*^9, 
   3.812744202921155*^9, 3.8127443024051905`*^9, 3.8127446461493416`*^9, 
   3.8127447325613327`*^9},
 CellLabel->
  "Out[300]=",ExpressionUUID->"5beaee48-0ef7-4136-9f4f-e7496fa19fd0"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Covariance", " ", "Calculation"}], ";", " ", 
    RowBox[{
     RowBox[{
      RowBox[{
      "this", " ", "is", " ", "one", " ", "calculation", " ", "approach"}], 
      "..."}], "another", " ", "is", " ", "below", " ", "and", " ", "yields", 
     " ", "the", " ", "same", " ", "result"}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"E", 
    RowBox[{"(", 
     RowBox[{"\[Pi]", "'"}], ")"}], "E", 
    RowBox[{"(", 
     RowBox[{"\[Pi]", "''"}], ")"}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"\[Pi]e1y", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"\[Pi]e1", "/.", 
        RowBox[{"y0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]e1y", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Pi]e1y", ",", "v", ",", "y0low", ",", "y0high"}], "]"}], 
      "*", "weighty"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"E\[Pi]e\[Theta]1", "=", 
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Pi]1\[Theta]", ",", "v", ",", "a1", ",", "b1"}], "]"}]}], 
     ";"}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]e2y", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"\[Pi]e2", "/.", 
        RowBox[{"y0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]e2y", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Pi]e2y", ",", "v", ",", "y0low", ",", "y0high"}], "]"}], 
      "*", "weighty"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"E\[Pi]e\[Theta]1", "=", 
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Pi]1\[Theta]", ",", "v", ",", "a1", ",", "b1"}], "]"}]}], 
     ";"}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]e12y", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"\[Pi]e1", "*", "\[Pi]e2"}], "/.", 
        RowBox[{"y0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]e12y", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Pi]e12y", ",", "v", ",", "y0low", ",", "y0high"}], "]"}], 
      "*", "weighty"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]e1\[Pi]", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Pi]e1y", "/.", 
        RowBox[{"\[Pi]0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]e1\[Pi]", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{
       "\[Pi]e1\[Pi]", ",", "v", ",", "\[Pi]0low", ",", "\[Pi]0high"}], "]"}],
       "*", "weight\[Pi]"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]e2\[Pi]", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Pi]e2y", "/.", 
        RowBox[{"\[Pi]0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]e2\[Pi]", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{
       "\[Pi]e2\[Pi]", ",", "v", ",", "\[Pi]0low", ",", "\[Pi]0high"}], "]"}],
       "*", "weight\[Pi]"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]e12\[Pi]", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Pi]e12y", "/.", 
        RowBox[{"\[Pi]0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]e12\[Pi]", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{
       "\[Pi]e12\[Pi]", ",", "v", ",", "\[Pi]0low", ",", "\[Pi]0high"}], 
       "]"}], "*", "weight\[Pi]"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]e1n", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Pi]e1\[Pi]", "/.", 
        RowBox[{"n0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]1e", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Pi]e1n", ",", "v", ",", "n0low", ",", "n0high"}], "]"}], 
      "*", "weightn"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]e2n", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Pi]e2\[Pi]", "/.", 
        RowBox[{"n0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]2e", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Pi]e2n", ",", "v", ",", "n0low", ",", "n0high"}], "]"}], 
      "*", "weightn"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Pi]e12n", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Pi]e12\[Pi]", "/.", 
        RowBox[{"n0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Pi]12e", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Pi]e12n", ",", "v", ",", "n0low", ",", "n0high"}], "]"}], 
      "*", "weightn"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{"covariance", " ", "=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"E\[Pi]12e", " ", "-", " ", 
       RowBox[{"E\[Pi]1e", "*", "E\[Pi]2e"}]}], ")"}], "/", "E\[Pi]2e", "  ", 
     RowBox[{"(*", 
      RowBox[{
      "this", " ", "is", " ", "covariance", " ", "of", " ", "alpha", " ", 
       "and", " ", "gamma"}], "*)"}]}]}]}]}]], "Input",
 CellChangeTimes->{{3.782142173537243*^9, 3.7821421821887627`*^9}, {
   3.7821423798666835`*^9, 3.7821423975959706`*^9}, {3.782143836016292*^9, 
   3.7821438866161613`*^9}, {3.7821439468314424`*^9, 3.78214415950089*^9}, {
   3.782144328814292*^9, 3.7821443486990566`*^9}, {3.7832558770327187`*^9, 
   3.783255880334742*^9}, {3.7832565644675684`*^9, 3.783256588129702*^9}, {
   3.7832566213652534`*^9, 3.7832566215051675`*^9}, {3.8049841943532033`*^9, 
   3.80498422209995*^9}, 3.8058958882876425`*^9, {3.8060648815734105`*^9, 
   3.806064897923929*^9}, {3.806157672762374*^9, 3.806157677438015*^9}, {
   3.8061609027009215`*^9, 3.8061609215951767`*^9}, {3.806161303849386*^9, 
   3.806161310345601*^9}, {3.8125891335203767`*^9, 3.812589153875243*^9}},
 CellLabel->
  "In[301]:=",ExpressionUUID->"e299b434-55cd-4bff-99ef-7c3fd0016d98"],

Cell[BoxData["0.5359616978986567`"], "Output",
 CellChangeTimes->{3.8058958906180887`*^9, 3.8060649677590966`*^9, 
  3.8061543171338725`*^9, 3.8061609337391386`*^9, 3.8061613134148254`*^9, 
  3.8127473945290613`*^9},
 CellLabel->
  "Out[319]=",ExpressionUUID->"b0fdeab6-1b29-40f8-bcb5-096be133b47e"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Components", " ", "of", " ", "\[CapitalDelta]"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"\[Alpha]", "=", 
     RowBox[{"\[Pi]e1", "-", "E\[Pi]1e"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Gamma]", "=", 
     RowBox[{"\[Pi]e2", "/", "E\[Pi]2e"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Sigma]2\[Alpha]y", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"\[Alpha]", "^", "2"}], "/.", 
        RowBox[{"y0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Sigma]2\[Alpha]y", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Sigma]2\[Alpha]y", ",", "v", ",", "y0low", ",", "y0high"}], 
       "]"}], "*", "weighty"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Sigma]2\[Alpha]\[Pi]", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Sigma]2\[Alpha]y", "/.", 
        RowBox[{"\[Pi]0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Sigma]2\[Alpha]\[Pi]", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{
       "\[Sigma]2\[Alpha]\[Pi]", ",", "v", ",", "\[Pi]0low", ",", 
        "\[Pi]0high"}], "]"}], "*", "weight\[Pi]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Sigma]2\[Alpha]n", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Sigma]2\[Alpha]\[Pi]", "/.", 
        RowBox[{"n0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Sigma]2\[Alpha]", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Sigma]2\[Alpha]n", ",", "v", ",", "n0low", ",", "n0high"}], 
       "]"}], "*", "weightn"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Sigma]\[Alpha]\[Gamma]y", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"\[Alpha]", "*", "\[Gamma]"}], "/.", 
        RowBox[{"y0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Sigma]\[Alpha]\[Gamma]y", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{
       "\[Sigma]\[Alpha]\[Gamma]y", ",", "v", ",", "y0low", ",", "y0high"}], 
       "]"}], "*", "weighty"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Sigma]\[Alpha]\[Gamma]\[Pi]", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Sigma]\[Alpha]\[Gamma]y", "/.", 
        RowBox[{"\[Pi]0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Sigma]\[Alpha]\[Gamma]\[Pi]", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{
       "\[Sigma]\[Alpha]\[Gamma]\[Pi]", ",", "v", ",", "\[Pi]0low", ",", 
        "\[Pi]0high"}], "]"}], "*", "weight\[Pi]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Sigma]\[Alpha]\[Gamma]n", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Sigma]\[Alpha]\[Gamma]\[Pi]", "/.", 
        RowBox[{"n0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Sigma]\[Alpha]\[Gamma]", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{
       "\[Sigma]\[Alpha]\[Gamma]n", ",", "v", ",", "n0low", ",", "n0high"}], 
       "]"}], "*", "weightn"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Gamma]2y", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"\[Gamma]", "^", "2"}], "/.", 
        RowBox[{"y0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Gamma]2y", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Gamma]2y", ",", "v", ",", "y0low", ",", "y0high"}], "]"}], 
      "*", "weighty"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Gamma]2\[Pi]", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Gamma]2y", "/.", 
        RowBox[{"\[Pi]0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Gamma]2\[Pi]", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{
       "\[Gamma]2\[Pi]", ",", "v", ",", "\[Pi]0low", ",", "\[Pi]0high"}], 
       "]"}], "*", "weight\[Pi]"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Gamma]2n", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Gamma]2\[Pi]", "/.", 
        RowBox[{"n0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Gamma]2", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Gamma]2n", ",", "v", ",", "n0low", ",", "n0high"}], "]"}], 
      "*", "weightn"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Sigma]2\[Gamma]", "=", 
     RowBox[{"E\[Gamma]2", "-", "1"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Alpha]\[Gamma]2y", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"\[Alpha]", "/", 
         RowBox[{"\[Gamma]", "^", "2"}]}], "/.", 
        RowBox[{"y0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Alpha]\[Gamma]2y", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Alpha]\[Gamma]2y", ",", "v", ",", "y0low", ",", "y0high"}], 
       "]"}], "*", "weighty"}]}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Alpha]\[Gamma]2\[Pi]", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Alpha]\[Gamma]2y", "/.", 
        RowBox[{"\[Pi]0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Alpha]\[Gamma]2\[Pi]", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{
       "\[Alpha]\[Gamma]2\[Pi]", ",", "v", ",", "\[Pi]0low", ",", 
        "\[Pi]0high"}], "]"}], "*", "weight\[Pi]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"\[Alpha]\[Gamma]2n", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"E\[Alpha]\[Gamma]2\[Pi]", "/.", 
        RowBox[{"n0", "\[Rule]", "#1"}]}], ")"}], "&"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"E\[Alpha]\[Gamma]2", "=", 
     RowBox[{
      RowBox[{"SimpsonsSum", "[", 
       RowBox[{"\[Alpha]\[Gamma]2n", ",", "v", ",", "n0low", ",", "n0high"}], 
       "]"}], "*", "weightn"}]}], ";"}], "\[IndentingNewLine]", 
   "\[IndentingNewLine]", "\"\<Covariance\>\"", "\[IndentingNewLine]", 
   "\[Sigma]\[Alpha]\[Gamma]", "\[IndentingNewLine]", "\[IndentingNewLine]", 
   "\"\<Adjustment to \[CapitalDelta]W\>\"", "\[IndentingNewLine]", 
   RowBox[{"W\[CapitalDelta]adjust", "=", 
    RowBox[{"\[Sigma]2\[Gamma]", "*", "\[CapitalDelta]W"}]}]}]}]], "Input",
 CellChangeTimes->{{3.8061543378847723`*^9, 3.80615435624647*^9}, {
  3.8061544210955734`*^9, 3.8061544353829412`*^9}, {3.806157932786916*^9, 
  3.806157972510036*^9}, {3.8061580817534313`*^9, 3.8061580954005456`*^9}, {
  3.806159038374971*^9, 3.8061593140112543`*^9}, {3.8061593517362847`*^9, 
  3.806159394893806*^9}, {3.8061598069565954`*^9, 3.80615981106621*^9}, {
  3.806159900380748*^9, 3.806159924890027*^9}, {3.8061602776921763`*^9, 
  3.806160281742626*^9}, {3.8061609722102284`*^9, 3.8061609988698435`*^9}, {
  3.8061610902037687`*^9, 3.8061612900113506`*^9}, {3.812589161710126*^9, 
  3.8125891656100736`*^9}},
 CellLabel->
  "In[320]:=",ExpressionUUID->"621c2a40-2155-4149-99a7-d3d8a1774cb4"],

Cell[BoxData["\<\"Covariance\"\>"], "Output",
 CellChangeTimes->{3.8061597746567793`*^9, 3.8061598309747295`*^9, 
  3.806159925714551*^9, 3.806160285161666*^9, 3.806161291254635*^9, 
  3.8127474019525833`*^9},
 CellLabel->
  "Out[347]=",ExpressionUUID->"84eed9ad-b45b-4dc7-9ab9-4bc2c15ee438"],

Cell[BoxData["0.5359616978986577`"], "Output",
 CellChangeTimes->{3.8061597746567793`*^9, 3.8061598309747295`*^9, 
  3.806159925714551*^9, 3.806160285161666*^9, 3.806161291254635*^9, 
  3.812747401955578*^9},
 CellLabel->
  "Out[348]=",ExpressionUUID->"5b946dc2-4255-4b2c-9311-48714993a582"],

Cell[BoxData["\<\"Adjustment to \[CapitalDelta]W\"\>"], "Output",
 CellChangeTimes->{3.8061597746567793`*^9, 3.8061598309747295`*^9, 
  3.806159925714551*^9, 3.806160285161666*^9, 3.806161291254635*^9, 
  3.812747401957575*^9},
 CellLabel->
  "Out[349]=",ExpressionUUID->"20fc62f6-6475-461c-9e5e-36dfb1d60855"],

Cell[BoxData["368470.63513354055`"], "Output",
 CellChangeTimes->{3.8061597746567793`*^9, 3.8061598309747295`*^9, 
  3.806159925714551*^9, 3.806160285161666*^9, 3.806161291254635*^9, 
  3.812747401959565*^9},
 CellLabel->
  "Out[350]=",ExpressionUUID->"19eaac78-c837-4221-bea3-401ab91e5e7b"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"\"\<Risk Adjustment\>\"", "\[IndentingNewLine]", 
   RowBox[{"Radjust", "=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "2"}], "*", 
       RowBox[{"d1", "/", 
        RowBox[{"(", 
         RowBox[{"E\[Pi]2e", "^", "2"}], ")"}]}]}], ")"}], "*", 
     "\[Sigma]2\[Alpha]", "*", "\[Sigma]2\[Gamma]"}]}], "\[IndentingNewLine]",
    "\[IndentingNewLine]", "\"\<Policy Adjustment\>\"", "\[IndentingNewLine]", 
   RowBox[{"Padjust", "=", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"\[Tau]1", "-", "\[Tau]1W"}], ")"}], "*", 
     RowBox[{"(", 
      RowBox[{"2", "*", 
       RowBox[{"d1", "/", 
        RowBox[{"(", 
         RowBox[{"2", "*", 
          RowBox[{"E\[Pi]2e", "^", "2"}]}], ")"}]}]}], ")"}], "*", 
     "E\[Alpha]\[Gamma]2"}]}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
   "\"\<Covariance Adjustment\>\"", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"\[CapitalDelta]true", "-", "\[CapitalDelta]W"}], ")"}], "-", 
    "W\[CapitalDelta]adjust", "-", "Radjust", "-", "Padjust"}], 
   "\[IndentingNewLine]"}]}]], "Input",
 CellChangeTimes->{{3.7821441346807933`*^9, 3.7821441355578966`*^9}, {
   3.8061594758648753`*^9, 3.806159535689961*^9}, {3.8061596119750423`*^9, 
   3.8061596908182135`*^9}, {3.806159727531517*^9, 3.8061597599033217`*^9}, {
   3.8061599601595645`*^9, 3.8061600292539597`*^9}, {3.806160327965672*^9, 
   3.806160331183941*^9}, {3.806160403555228*^9, 3.8061604059740534`*^9}, {
   3.8061614325444856`*^9, 3.8061614336678467`*^9}, {3.806161463686601*^9, 
   3.806161464856926*^9}, 3.8061616049946785`*^9},
 CellLabel->
  "In[351]:=",ExpressionUUID->"b207a403-e7ad-4a38-b6b6-288380782685"],

Cell[BoxData["\<\"Risk Adjustment\"\>"], "Output",
 CellChangeTimes->{
  3.8061598426640053`*^9, 3.806159928920703*^9, {3.806160024881463*^9, 
   3.8061600300015326`*^9}, 3.806160428625193*^9, 3.8061614669257936`*^9, 
   3.806161634450798*^9, 3.812747407646471*^9},
 CellLabel->
  "Out[351]=",ExpressionUUID->"b2207a3e-186c-4b97-982b-1b763dda3694"],

Cell[BoxData[
 RowBox[{"-", "1.6754554923809033`*^6"}]], "Output",
 CellChangeTimes->{
  3.8061598426640053`*^9, 3.806159928920703*^9, {3.806160024881463*^9, 
   3.8061600300015326`*^9}, 3.806160428625193*^9, 3.8061614669257936`*^9, 
   3.806161634450798*^9, 3.8127474076544743`*^9},
 CellLabel->
  "Out[352]=",ExpressionUUID->"ecf7dde0-e37b-4dc5-b010-25850f733b6e"],

Cell[BoxData["\<\"Policy Adjustment\"\>"], "Output",
 CellChangeTimes->{
  3.8061598426640053`*^9, 3.806159928920703*^9, {3.806160024881463*^9, 
   3.8061600300015326`*^9}, 3.806160428625193*^9, 3.8061614669257936`*^9, 
   3.806161634450798*^9, 3.8127474076564674`*^9},
 CellLabel->
  "Out[353]=",ExpressionUUID->"ec54d87c-74b7-4d69-8385-c28332fd2b1a"],

Cell[BoxData["1.0492183550677825`*^7"], "Output",
 CellChangeTimes->{
  3.8061598426640053`*^9, 3.806159928920703*^9, {3.806160024881463*^9, 
   3.8061600300015326`*^9}, 3.806160428625193*^9, 3.8061614669257936`*^9, 
   3.806161634450798*^9, 3.8127474076584377`*^9},
 CellLabel->
  "Out[354]=",ExpressionUUID->"eeceaf41-1011-4b35-a937-5826fc3804f1"],

Cell[BoxData["\<\"Covariance Adjustment\"\>"], "Output",
 CellChangeTimes->{
  3.8061598426640053`*^9, 3.806159928920703*^9, {3.806160024881463*^9, 
   3.8061600300015326`*^9}, 3.806160428625193*^9, 3.8061614669257936`*^9, 
   3.806161634450798*^9, 3.812747407660467*^9},
 CellLabel->
  "Out[355]=",ExpressionUUID->"6e0e10ef-a446-4fbd-b694-b5993cbe7c8f"],

Cell[BoxData[
 RowBox[{"-", "1.3414633333422042`*^7"}]], "Output",
 CellChangeTimes->{
  3.8061598426640053`*^9, 3.806159928920703*^9, {3.806160024881463*^9, 
   3.8061600300015326`*^9}, 3.806160428625193*^9, 3.8061614669257936`*^9, 
   3.806161634450798*^9, 3.8127474076634245`*^9},
 CellLabel->
  "Out[356]=",ExpressionUUID->"28b8943c-af99-4010-a25f-736fa897beae"]
}, Open  ]]
},
WindowSize->{1280, 941},
WindowMargins->{{-8, Automatic}, {Automatic, -8}},
Magnification->1.25,
FrontEndVersion->"12.0 for Microsoft Windows (64-bit) (April 8, 2019)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 883, 20, 108, "Input",ExpressionUUID->"851da910-5bec-44c6-af44-3b60a1f63ce3"],
Cell[1444, 42, 10483, 203, 612, "Input",ExpressionUUID->"395c95fd-1bbf-41d2-8e50-37061b2a08ce"],
Cell[CellGroupData[{
Cell[11952, 249, 2730, 64, 204, "Input",ExpressionUUID->"daf9159f-7516-4d32-b91a-decf05d4fe44"],
Cell[14685, 315, 1262, 22, 29, "Message",ExpressionUUID->"b8ac0353-cf74-49a5-b1ed-89d483392b61"]
}, Open  ]],
Cell[15962, 340, 11058, 200, 588, "Input",ExpressionUUID->"7b96dfab-41a9-4b2b-a3f4-5f1c44b9c1db"],
Cell[27023, 542, 10491, 177, 693, "Input",ExpressionUUID->"71d7e32a-4ae8-4ebb-b902-48a504476939"],
Cell[37517, 721, 8499, 191, 756, "Input",ExpressionUUID->"cbf0001c-ea66-4e0d-ac41-d76bb60a4987"],
Cell[46019, 914, 5159, 118, 468, "Input",ExpressionUUID->"1a7c9eee-15a8-4798-a586-1912b00536eb"],
Cell[51181, 1034, 4961, 111, 444, "Input",ExpressionUUID->"23f8ba2f-2397-4804-9220-95de22826350"],
Cell[56145, 1147, 6288, 158, 852, "Input",ExpressionUUID->"8f646553-56c0-4a66-9846-7ab75b07e623"],
Cell[62436, 1307, 7174, 186, 900, "Input",ExpressionUUID->"ea5b2a10-3df4-4037-83d4-434ad2e8e09a"],
Cell[CellGroupData[{
Cell[69635, 1497, 8253, 203, 1692, "Input",ExpressionUUID->"9a82675f-d59c-43ec-99e2-4b1b09762af4"],
Cell[77891, 1702, 760, 12, 40, "Output",ExpressionUUID->"7423f4b1-0b2c-4169-82f0-4ec57d895ff3"],
Cell[78654, 1716, 757, 12, 40, "Output",ExpressionUUID->"0c6a1c50-313e-419b-a82c-58af9a48579b"],
Cell[79414, 1730, 757, 12, 40, "Output",ExpressionUUID->"16b0601f-0310-4a4e-975d-440d4537ffc9"],
Cell[80174, 1744, 745, 12, 40, "Output",ExpressionUUID->"d2b04fa4-20af-46c4-bfe0-005335eee82f"],
Cell[80922, 1758, 799, 13, 40, "Output",ExpressionUUID->"cd19888a-8180-4c46-aa0b-b6cde8cfb859"],
Cell[81724, 1773, 753, 12, 40, "Output",ExpressionUUID->"42092a27-8fb1-4b10-8c63-e3ed52c27d87"],
Cell[82480, 1787, 745, 12, 40, "Output",ExpressionUUID->"4c44ea21-4065-4b39-b0fd-07f64174236d"],
Cell[83228, 1801, 762, 12, 40, "Output",ExpressionUUID->"c0f13d81-c851-4863-9857-8f2957d2c7c1"],
Cell[83993, 1815, 755, 12, 40, "Output",ExpressionUUID->"9f4dc2a3-9d73-45f8-a786-bf7f2c728f2c"],
Cell[84751, 1829, 743, 12, 40, "Output",ExpressionUUID->"2df40f5c-4a1d-4dc0-9d58-d88a62e987c8"],
Cell[85497, 1843, 801, 13, 40, "Output",ExpressionUUID->"ce566f49-dc23-40fc-9ccf-7d417c0f2a92"],
Cell[86301, 1858, 756, 12, 40, "Output",ExpressionUUID->"2c91a6a1-ba51-490f-abb9-433eaa74e5d9"],
Cell[87060, 1872, 743, 12, 40, "Output",ExpressionUUID->"10be7a1f-c009-4b1e-8753-a96232f0a8bd"],
Cell[87806, 1886, 757, 12, 40, "Output",ExpressionUUID->"b674b698-2da9-4dba-89b0-a34c9b946fcf"],
Cell[88566, 1900, 754, 12, 40, "Output",ExpressionUUID->"6a5a2ffe-bd7b-47db-a82f-b5b7b5a2ccf1"],
Cell[89323, 1914, 743, 12, 40, "Output",ExpressionUUID->"d6061b84-bfab-4d2e-972a-1c85c5a33594"],
Cell[90069, 1928, 765, 12, 40, "Output",ExpressionUUID->"84c0e03c-22cb-43f3-84ff-c16d735063d1"],
Cell[90837, 1942, 772, 13, 40, "Output",ExpressionUUID->"3db2f32c-5c86-4bcc-a592-538263016639"],
Cell[91612, 1957, 743, 12, 40, "Output",ExpressionUUID->"98b3ecb8-a03e-4c41-96d7-1fe52fdfa3e1"],
Cell[92358, 1971, 775, 12, 40, "Output",ExpressionUUID->"293c2cb0-fcb0-4dce-b369-406e4ea70790"],
Cell[93136, 1985, 745, 12, 40, "Output",ExpressionUUID->"22544025-0796-45e2-a9a4-45b7bf9ce576"],
Cell[93884, 1999, 745, 12, 40, "Output",ExpressionUUID->"33ed285a-78bb-4ce2-9826-fac27fdac0ef"],
Cell[94632, 2013, 761, 12, 40, "Output",ExpressionUUID->"1aeaa1b7-e00d-47db-807b-01fb29d11da0"],
Cell[95396, 2027, 755, 12, 40, "Output",ExpressionUUID->"39d8c8e2-1bc4-4d44-b902-ac9bbe1ed20b"],
Cell[96154, 2041, 743, 12, 40, "Output",ExpressionUUID->"4a3ac2f6-c934-4bf2-b198-f9c4f6acede0"],
Cell[96900, 2055, 823, 13, 40, "Output",ExpressionUUID->"5fc91c3f-8abd-477d-ba78-4192d35e7eb2"],
Cell[97726, 2070, 769, 13, 40, "Output",ExpressionUUID->"31d895ef-e65d-47bf-a50d-53dcb76d6036"],
Cell[98498, 2085, 745, 12, 40, "Output",ExpressionUUID->"2a4ead87-757e-4e7e-9ca0-c893d558f881"],
Cell[99246, 2099, 762, 12, 40, "Output",ExpressionUUID->"598fbaf1-4ec0-4008-ac8f-c1394206183e"],
Cell[100011, 2113, 758, 12, 40, "Output",ExpressionUUID->"bf5aea58-8446-4715-89c5-a0d1e5989cbd"],
Cell[100772, 2127, 745, 12, 40, "Output",ExpressionUUID->"d4e1a51b-8b30-40fb-b473-a98c80a0e592"],
Cell[101520, 2141, 801, 13, 40, "Output",ExpressionUUID->"c9ffc715-cc93-446d-8c6f-45a67a00d250"],
Cell[102324, 2156, 757, 12, 40, "Output",ExpressionUUID->"a81caffd-af12-4f8d-a8cc-61acfb99eaef"],
Cell[103084, 2170, 742, 12, 40, "Output",ExpressionUUID->"0978c7b0-870c-44da-94e6-354033661617"],
Cell[103829, 2184, 764, 12, 40, "Output",ExpressionUUID->"b7684e05-7437-4101-89b0-013867be5357"],
Cell[104596, 2198, 773, 13, 40, "Output",ExpressionUUID->"88fe9d5d-c381-459e-8f50-977b402a0920"],
Cell[105372, 2213, 743, 12, 40, "Output",ExpressionUUID->"e13a8f71-755b-47e2-9a6a-5cbd3c90dcec"],
Cell[106118, 2227, 745, 12, 40, "Output",ExpressionUUID->"bddb10ef-6eba-4142-bdd1-161741289484"],
Cell[106866, 2241, 805, 13, 40, "Output",ExpressionUUID->"34a2ee8b-d8b5-4483-9ba2-7b37830629fd"],
Cell[107674, 2256, 755, 12, 40, "Output",ExpressionUUID->"13a7ab4e-e840-4c2e-a6db-44f650e1add2"],
Cell[108432, 2270, 745, 12, 40, "Output",ExpressionUUID->"22679df2-3f69-4121-aad1-16fa555915d1"],
Cell[109180, 2284, 764, 12, 40, "Output",ExpressionUUID->"995b8406-4a7b-4f4d-b0a1-1a14e364696e"],
Cell[109947, 2298, 751, 12, 40, "Output",ExpressionUUID->"4514676c-be1f-43dd-ab07-0b8fdf8cb745"],
Cell[110701, 2312, 745, 12, 40, "Output",ExpressionUUID->"b5557ed8-8252-4ff9-8f82-52f3cd06ed5e"],
Cell[111449, 2326, 752, 12, 40, "Output",ExpressionUUID->"811e8c7e-2750-4d39-b0a0-9030ef15523f"],
Cell[112204, 2340, 757, 12, 40, "Output",ExpressionUUID->"c6a5603d-75e4-4322-84be-395c6d1ebf07"],
Cell[112964, 2354, 743, 12, 40, "Output",ExpressionUUID->"189ab2bc-102e-4d56-b4ac-0541b23883cf"],
Cell[113710, 2368, 858, 13, 40, "Output",ExpressionUUID->"2f9759df-b04a-4072-8cb4-6c83fb4eb671"],
Cell[114571, 2383, 758, 12, 40, "Output",ExpressionUUID->"fe0bd304-f75e-4089-9d18-9ae8266cf8da"],
Cell[115332, 2397, 743, 12, 40, "Output",ExpressionUUID->"1751448f-97f4-4bbe-9ce6-7ff470dd96ff"],
Cell[116078, 2411, 763, 12, 40, "Output",ExpressionUUID->"a6a27a59-05a5-4d94-a8b7-6c5602a41c41"],
Cell[116844, 2425, 756, 12, 40, "Output",ExpressionUUID->"8ac3ab23-c685-49dc-ba11-48859476cc02"],
Cell[117603, 2439, 745, 12, 40, "Output",ExpressionUUID->"7f40b1dd-b8d2-4a48-becb-25c77a04e709"],
Cell[118351, 2453, 745, 12, 40, "Output",ExpressionUUID->"244f61c9-f2a9-4cf6-835d-0b3d63ed0f72"],
Cell[119099, 2467, 792, 13, 40, "Output",ExpressionUUID->"ac00df29-b43d-4af8-9791-7905764d6173"],
Cell[119894, 2482, 755, 12, 40, "Output",ExpressionUUID->"30354b1b-8066-4d60-9e61-75232f5736cb"],
Cell[120652, 2496, 743, 12, 40, "Output",ExpressionUUID->"06aefbb0-befe-4091-b03e-4afbd4579a93"],
Cell[121398, 2510, 765, 12, 40, "Output",ExpressionUUID->"a0c3ca86-63c9-4ef5-8159-e209f3d1f76c"],
Cell[122166, 2524, 775, 13, 40, "Output",ExpressionUUID->"2bd99fcb-b5a5-4aab-b064-7818f76cc0d4"],
Cell[122944, 2539, 745, 12, 40, "Output",ExpressionUUID->"dec9bfbf-399e-4e94-b63d-10f41c8dfad2"],
Cell[123692, 2553, 808, 13, 40, "Output",ExpressionUUID->"e120990f-d8ff-46b7-bc5c-7720d87511f6"],
Cell[124503, 2568, 758, 12, 40, "Output",ExpressionUUID->"8f84a9cf-4759-4a27-a3c2-c44650725801"],
Cell[125264, 2582, 745, 12, 40, "Output",ExpressionUUID->"8582eecf-4b95-4f34-b44e-38f3c865f194"],
Cell[126012, 2596, 804, 13, 40, "Output",ExpressionUUID->"d563a08f-1f0e-4455-b259-652a071f8b2a"],
Cell[126819, 2611, 753, 12, 40, "Output",ExpressionUUID->"82d910d7-27ed-4562-8fda-ce6634e2447e"],
Cell[127575, 2625, 810, 13, 40, "Output",ExpressionUUID->"7535d366-2349-41ae-8e6e-00c9b4b807a3"],
Cell[128388, 2640, 755, 12, 40, "Output",ExpressionUUID->"5beaee48-0ef7-4136-9f4f-e7496fa19fd0"]
}, Open  ]],
Cell[CellGroupData[{
Cell[129180, 2657, 6477, 175, 780, "Input",ExpressionUUID->"e299b434-55cd-4bff-99ef-7c3fd0016d98"],
Cell[135660, 2834, 299, 5, 40, "Output",ExpressionUUID->"b0fdeab6-1b29-40f8-bcb5-096be133b47e"]
}, Open  ]],
Cell[CellGroupData[{
Cell[135996, 2844, 7621, 204, 948, "Input",ExpressionUUID->"621c2a40-2155-4149-99a7-d3d8a1774cb4"],
Cell[143620, 3050, 292, 5, 40, "Output",ExpressionUUID->"84eed9ad-b45b-4dc7-9ab9-4bc2c15ee438"],
Cell[143915, 3057, 291, 5, 40, "Output",ExpressionUUID->"5b946dc2-4255-4b2c-9311-48714993a582"],
Cell[144209, 3064, 310, 5, 40, "Output",ExpressionUUID->"20fc62f6-6475-461c-9e5e-36dfb1d60855"],
Cell[144522, 3071, 291, 5, 40, "Output",ExpressionUUID->"19eaac78-c837-4221-bea3-401ab91e5e7b"]
}, Open  ]],
Cell[CellGroupData[{
Cell[144850, 3081, 1730, 38, 252, "Input",ExpressionUUID->"b207a403-e7ad-4a38-b6b6-288380782685"],
Cell[146583, 3121, 348, 6, 40, "Output",ExpressionUUID->"b2207a3e-186c-4b97-982b-1b763dda3694"],
Cell[146934, 3129, 366, 7, 40, "Output",ExpressionUUID->"ecf7dde0-e37b-4dc5-b010-25850f733b6e"],
Cell[147303, 3138, 352, 6, 40, "Output",ExpressionUUID->"ec54d87c-74b7-4d69-8385-c28332fd2b1a"],
Cell[147658, 3146, 349, 6, 40, "Output",ExpressionUUID->"eeceaf41-1011-4b35-a937-5826fc3804f1"],
Cell[148010, 3154, 354, 6, 40, "Output",ExpressionUUID->"6e0e10ef-a446-4fbd-b694-b5993cbe7c8f"],
Cell[148367, 3162, 366, 7, 73, "Output",ExpressionUUID->"28b8943c-af99-4010-a25f-736fa897beae"]
}, Open  ]]
}
]
*)

